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REAL-AIR DATA-REDUC TION PROCEDURES 
BASED ON FLOW PARAMETERS MEASURED IN THE TEST SECTION 
OF SUPERSONIC AND HYPERSONIC FACILITIES 

By Charles G. Miller HI and Sue E. Wilder 
Langley Research Center 

SUMMARY 

Data- reduction procedures for determining free-stream and post-normal-shock 
kinetic and thermodynamic quantities are derived. These procedures are applicable to 
imperfect real-air flows in thermochemical equilibrium for temperatures to 15 000 K 
and a range of pressures from 0.25 N/m2 to 1 GN/m2. Although derived primarily to 
meet the immediate needs of the Langley 6 -inch expansion tube, these procedures are 
applicable to any supersonic or hypersonic real-air test facility where combinations of 
the following flow parameters are measured in the test section: 

(1) Stagnation pressure behind normal shock 

(2) Free-stream static pressure 

(3) Stagnation-point heat-transfer rate 

(4) Free-stream flow velocity 

(5) Stagnation density behind normal shock 

(6) Free-stream density 

Nine data- reduction procedures resulting from various combinations of three of these 
measured flow parameters are derived. These procedures employ an adjustment of 
computed flow parameters by numerical iteration until measured and computed flow 
parameters are within a prescribed tolerance. 

Because the above six flow parameters are measured in the test section, these 
procedures do not depend explicitly upon measured or calculated upstream flow param- 
eters. The elimination of dependence on upstream flow conditions results in a reduction 
in the uncertainty in predicted test-section conditions. 

Limitations of the various procedures and uncertainties in calculated flow quantities 
corresponding to uncertainties in measured input data are discussed. All nine procedures 
are incorporated into a single computer program written in FORTRAN IV language. A 
listing of this computer program is presented, along with a description of the inputs 
required and a sample of the data printout. 



INTRODUCTION 


Over the past decade, a number of studies (refs. 1 to 11) have been directed toward 
prediction of performance characteristics of expansion tubes and expansion tunnels. The 
initial theoretical study of the expansion tube was performed by Trimpi (ref. 1). In refer- 
ence 1 a simplified flow model based on idealized processes was used. However, Trimpi 
acknowledged the possible existence of detrimental effects on expansion tube performance 
arising from real physical conditions (such as noninstantaneous primary- and secondary- 
diaphragm rupture, shock-wave attenuation, interface mixing, and flow turbulence). At 
the time reference 1 was written, the extent to which these possible phenomena deviated 
the actual flow from idealized flow remained to be determined experimentally. (The 
complexity associated with these real physical conditions generally prohibits rigorous 
theoretical treatment.) To this end, the Langley shock tunnel was modified in 1961 to 
serve as a pilot model expansion tube. Experimental results of exploratory studies in 
this facility are reported in references 12 and 13. As might be expected, significant dif- 
ferences between measured and theoretically predicted flow quantities were observed. 

Investigations representing extensions of existing experimental studies into the 
high-velocity real-air regime will be performed in the Langley 6 -inch expansion tube. 

To make meaningful comparisons of these real-air data with the existing data and corre- 
sponding correlations, accurate predictions of free-stream and post-normal-shock flow 
conditions are required. 

The usual procedure for determining expansion tube test-section flow conditions 
(for real air in thermochemical equilibrium) is first to determine conditions behind the 
incident shock propagating into the static test gas within the intermediate section (see 
fig. 1). (The charts of ref. 14 are a convenient and often employed means of determining 
these conditions.) An isentropic unsteady expansion from the velocity of the shocked test 
gas in the intermediate section to the measured free-stream velocity at the test section 
is performed (ref. 3). With the free-stream flow conditions determined, the post-normal- 
shock conditions are determined from existing tables or charts for standing normal shocks. 
This scheme generally yields test- section flow quantities significantly different from mea- 
sured quantities. 

The purpose of this study is to provide a means for obtaining accurate test- section 
flow conditions in the expansion tube and expansion tunnel. The computational scheme 
used to predict test-section conditions is based on flow parameters measured in the test 
section, and thus an explicit dependence upon measured or calculated upstream flow 
parameters is eliminated. The elimination of dependence on upstream flow conditions 
should result in a substantial reduction in the uncertainty in predicted test-section condi- 
tions. (For example, it has been speculated in refs. 13 and 15 that one cause of failure 
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of the theoretical approach in predicting test-section conditions is nonideal rupture of 
the upstream secondary diaphragm.) Such a computational scheme requires measure- 
ment of three flow quantities at the test section in order to satisfy the conservation equa- 
tions for a standing normal shock. Presently, three test-section flow quantities are mea- 
sured, on a routine basis, in the Langley 6 -inch expansion tube: 

(1) Stagnation pressure behind normal shock 

(2) Tube wall static pressure 

(3) Velocity of interface of acceleration gas and test gas 

Another measurable flow quantity that will be obtained on a routine basis is stagnation- 
point heat-transfer rate. Stagnation-point heat-transfer rate measurements have been 
obtained successfully in shock tubes and shock tunnels with thin-film gages (e.g., refs. 16 
to 18) and thick-film gages (e.g., refs. 19 and 20). Hence, these instrumentation tech- 
niques are applicable in the expansion tube (ref. 21). Stagnation-point heating rate has 
been employed as a basic input datum in data-reduction procedures for arc-heated 
impulse tunnels (commonly referred to as hotshot tunnels) as discussed in reference 22. 
Following the example of reference 22, the stagnation-point heating rate is included 
herein as a fourth basic input datum. 

An additional flow quantity that can be inferred from experimental techniques is the 
density. Although free-stream density and stagnation -point density behind a normal 
shock are not presently measured in the Langley 6 -inch expansion tube or expansion tun- 
nel on a routine basis, these two quantities are, nevertheless, considered herein as input 
data. 

From these six flow quantities measured at the test section, nine procedures for 
determining free-stream and post-normal-shock flow conditions are derived. These 
procedures utilize the thermodynamic properties for imperfect real air in thermochemi- 
cal equilibrium as tabulated in references 23 and 24. (The results of refs. 23 and 24 
represent a compilation of the air data of refs. 25 and 26, with interpolations and differ- 
entiation. They cover a range of temperatures from 100 K to 15 000 K and a range of 
pressures from 0.25 N/m^ to 1 GN/m2.) Although derived to meet the immediate needs 
of the Langley 6-inch expansion tube, these procedures are also applicable to supersonic 
or hypersonic test facilities with real-air flows. 

Limitations of the procedures are discussed. An experimental uncertainty is 
assigned to each of the six experimental inputs and the corresponding uncertainty in cal- 
culated flow quantities is examined. 

The relations considered for the prediction of stagnation-point heat-transfer rate 
for real air are discussed in appendix A. The transport properties for high-temperature 
air required in the prediction of stagnation-point heat-transfer rate are discussed in 
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appendix B. The calculation schemes for the individual data- reduction procedures are 
presented in appendix C. All nine procedures are incorporated into a single computer 
program written in FORTRAN IV language. A listing of this computer program and a 
sample data printout are presented in appendix D, and a description of the inputs required 
is presented in appendix E. 


* SYMBOLS 

The International System of Units (SI) is used for all physical quantities in this 
study. Conversion factors relating the SI Units to U.S. Customary Units are given in 
reference 27. 

a speed of sound, m/sec 

Cp specific heat at constant pressure, J/kg-K 

D diffusion coefficient, m^/sec 

h specific enthalpy, m2/sec2 (j/kg) 

he convective heat- transfer coefficient, W/m2-K 

h n atomic dissociation energy, J/kg 


k thermal conductivity, W/m-K 

M Mach number, V/a 

Lewis number, pc^o/k 
n Nu Nusselt number, h^s/k 

Np r Prandtl number, pc p y/k 

N Re 
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Reynolds number per meter, pV/p. 
pressure, N/m2 
heat-transfer rate, W/m2 



universal gas constant, 8.31434 kJ/kmol-K 


R 

r e 

r g 

S 

s/R 

T 

V 
W 

W 0 

z 

z* 

/3 

V s PCpk 


y E 


P 


P 

T 


effective nose radius, m 
geometric nose radius, m 

distance along body measured from stagnation point, m 

nondimensional specific entropy 

temperature, K 

velocity, m/sec 

molecular weight, kg/kmol 

molecular weight of undissociated air, 28.967 kg/kmol 
compressibility factor, pW 0 /^RT 

ratio of number of moles to number of moles for undissociated air, W 0 /w 
stagnation- point velocity gradient, sec - l 


isentropic exponent, 

V 9 lQ g PJs/R 

coefficient of viscosity, N-sec/m2 

nondimensionalized ratio of uncertainty in calculated flow quantity to corre- 
sponding uncertainty in input quantity (see eq. (9)) 

density, kg/m3 

time, sec 
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Subscripts: 


amb ambient 

av average 

c calculated 

low lower limit 

m measured 

prev previous value of a parameter 

5 based on distance along body from stagnation point 

s model surface material 

t stagnation conditions behind normal shock 

up upper limit 

w model wall 

w,l tube or nozzle wall 

1 free stream 

2 static conditions immediately behind normal shock 

Superscripts: 

a order of iteration 

~ approximate value 

ANALYSIS 

Several topics are presented before the discussion of the various calculation pro- 
cedures for determining free-stream and post-normal-shock flow quantities. First a 
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brief description of the Langley 6-inch expansion tube and expansion tunnel is given. This 
description is followed by a discussion of the instrumentation techniques employed to 
obtain the experimental data input required by the various calculation procedures. The 
magnitude of the uncertainties in the experimental data input is also discussed. The 
source of the thermodynamic properties for real air in thermochemical equilibrium 
employed in this study is discussed. Next, the methods for crossing the standing normal 
shock are discussed, since these methods are common to all calculation procedures. 
Finally, the relations used to predict stagnation-point heat-transfer rate are discussed 
briefly. 


Description of Expansion Tube and Expansion Tunnel 

The Langley 6-inch expansion tube is basically a cylindrical tube with a 15.24-cm 
inside diameter, divided by two diaphragms into three sections. The most upstream sec- 
tion is the driver or high-pressure section. This section is pressurized at ambient tem- 
perature with a gas having a high speed of sound, such as hydrogen or helium. (Greater 
operation efficiency is realized with gases having a high speed of sound; e.g., see ref. 28.) 
The pressure and speed of sound of the driver gas are increased further by heating the 
gas with a 3 -MW resistance heater or utilizing an arc discharge into the gas from a 10-MJ 
capacitor bank. The intermediate section is usually referred to as the driven section. 

This section is evacuated and filled with the test gas at ambient temperature. The most 
downstream section is referred to as the expansion or acceleration section. This section 
is also evacuated and generally filled with helium at a low pressure and ambient tempera- 
ture. For resistance heating of the driver gas, the driver and driven sections are sepa- 
rated by a double diaphragm apparatus capable of withstanding a maximum pressure dif- 
ferential of 68.95 MN/m2. (By controlling the pressure level in the small chamber 
between these diaphragms, the time of diaphragm rupture can be controlled.) For arc 
heating, a single diaphragm is used between the driver and driven sections. A weak, low- 
pressure diaphragm (secondary diaphragm) separates the driven and acceleration sections. 
The test section and model are located at the downstream exit of the acceleration section. 

The operating sequence, which is shown schematically in figure 1, begins with the 
rupture of the high-pressure diaphragm. A primary shock wave propagates into the static 
test gas and an expansion wave propagates into the driver gas. The shock wave then 
encounters and ruptures the low-pressure diaphragm. A secondary shock wave propagates 
into the low-pressure accelerating gas while an upstream expansion wave moves into the 
test gas. In passing through this upstream expansion wave, which is being washed down- 
stream since the shock-heated test gas is supersonic, the test gas undergoes an isentropic 
unsteady expansion resulting in an increase in the flow velocity. 

The expansion tunnel is simply an expansion tube with a nozzle added at the down- 
stream end. Thus, the test gas is processed first by the primary shock in the interme- 
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diate section, then by an unsteady expansion in the acceleration chamber, and finally by 
an isentropic steady expansion in the nozzle (ref. 4). 

Experimental Data 

Pressure .- The expansion tube is characterized by extremely short test times. In 
general, the test time is less than 400 psec, and thus the pressure instrumentation must 
have very fast response to pressure change and a minimum of orifice -cavity volume to 
reduce pressure-lag effects. Presently, stagnation pressure behind a normal shock (pitot 
pressure) and expansion tube or nozzle-wall pressure are measured with miniature piezo- 
electric (quartz) transducers having rise times of approximately 1 to 3 [isec and a pres- 
sure range of approximately 700 N/m2 to 20 MN/m2. The pressure transducers are used 
in conjunction with a charge amplifier, and the output signal is recorded from an oscillo- 
scope with the aid of a camera. 

Experimental uncertainties in such pressure measurements are dependent on many 
factors, such as pressure level with regard to transducer sensitivity, calibration technique 
(static or dynamic), transducer linearity, oscilloscope accuracy, quality of oscilloscope 
traces with respect to the signal-to-noise ratio, and data read-up procedure. Hence, a 
general assignment of the uncertainty in these pressure measurements is not possible. 

On the basis of previous experience, the maximum uncertainties in pressure measure- 
ments are believed to be less than ±20 percent for tube or nozzle-wall pressure and for 
pitot pressure. 

Velocity .- The free-stream velocity is usually inferred by using a microwave tech- 
nique. A microwave signal is propagated upstream into the acceleration section by means 
of an antenna mounted in the neighborhood of the acceleration-section exit (test section). 
The microwave signal is reflected from the interface, provided the electron concentration 
at the interface is high enough for reflection of the operating signal. (As shown in ref. 9, 
for incident shock velocities less than 10 km/sec in the helium acceleration gas, the 
helium is not ionized and is therefore transparent to the microwave signal.) The inter- 
ference of the reflected wave with the transmitted wave provides a measurement of the 
interface velocity. The microwave signal is recorded on film by means of a combination 
of a high-speed drum camera and an oscilloscope. The uncertainty in measuring the 
interface velocity with this technique is believed to be less than ±2 percent. However, 
the experimental results of reference 29 show that the flow velocity is not uniform behind 
the interface, and inferring that the free-stream velocity is equal to the interface velocity 
may result in errors up to 5 percent for velocities around 6 km/sec. For the purposes 
of this study, the uncertainty in free-stream velocity (as inferred from microwave mea- 
surements) will be assumed not to exceed ±5 percent. 
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Stagnation- point heat- transfer rate .- Because of the heat-flux limitations of thin- 
film heat-transfer gages, thick- film or calorimeter heat-transfer gages are used to mea- 
sure stagnation-point heat-transfer rates in the Langley 6-inch expansion tube. The thick- 
film gage consists of a thin-foil sensing element (usually pure platinum) mounted on an 
insulating substrate. The foil thickness is chosen so that only a negligible amount of heat 
is transferred to the substrate (generally less than 5 percent) during the testing period. 

A detailed theoretical treatment of the thick-film heat-transfer gage and data- reduction 
procedure is presented in reference 19. In reference 20, the experimental uncertainty 
for the thick-film heat-transfer gage is stated to be ±11 percent. Possible detrimental 
effects due to flow turbulence level and flow contamination, the magnitudes of which are 
unknown for the Langley 6 -inch expansion tube, would tend to increase the uncertainty in 
stagnation- point heat-transfer measurements. For this study, the uncertainty in these 
measurements will be assumed not to exceed ±20 percent. 

Density . - Free-stream density measurements have been made in the Langley pilot 
model expansion tube by using a spectroscopic method similar to that described in refer- 
ence 30. These measurements, which are based on the light- absorption properties of 
oxygen molecules, are believed to have an uncertainty of approximately ±10 percent. For 
purposes of this study, the uncertainty in both free-stream density and stagnation-point 
density will be assumed not to exceed ±20 percent. 

Thermodynamic Properties for Real Air 

Thermodynamic properties for imperfect real air in thermochemical equilibrium 
are obtained from a magnetic tape furnished to the Langley Research Center by the Arnold 
Engineering Development Center (AEDC) in late 1965. The thermodynamic properties 
obtained from this tape correspond to the properties tabulated in reference 23 for various 
values of s/R. The temperature range of reference 23 and of the AEDC tape is 100 K to 
15 000 K and the pressure range is 0.25 N/m2 to 1 GN/m2. For temperatures from 100 K 
to 1500 K, the results of reference 23 were obtained by interpolation of the data of refer- 
ence 25; for temperatures from 1500 K to 15 000 K, the results were obtained by interpo- 
lation of the data of reference 26. In addition to the effects of dissociation and ionization, 
the high- temperature data of reference 26 include second virial corrections for interac- 
tions between neutral-neutral species and ion-neutral species. The low-temperature 
(undissociated) air composition of reference 26 is 78.084 percent N2, 20.946 percent O2, 
0.934 percent Ar, 0.033 percent CO2, and 0.003 percent Ne by volume. This composition 
corresponds to a molecular weight of 28.967. 

Thermodynamic properties included on the composite AEDC tape are a, h, p, 
s/R, T, Z, Z* , y^, and p. The properties a and y^, correspond to those tabulated 
in reference 31 for temperatures from 1500 K to 15 000 K. The results of reference 31 
were obtained by interpolation and differentiation of the real-air data of reference 26. 
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As obtained from AEDC, the subroutine for searching the real-air tape required 
inputs of s/R in conjunction with one of the following thermodynamic properties: a, h, 
p, T, p. (The procedure whereby a, h, p, T, or p is used in conjunction with s/R 
for inputs to the AEDC tape is designated herein as SLOW.) An interpolation procedure 
allowing combinations of h, p, and p as inputs to the AEDC tape was derived for this 
study. This interpolation procedure is referred to herein as SEARCH (L), where L = 1 
denotes inputs h and p, L = 2 denotes inputs p and p, and L = 3 denotes inputs 
h and p. 

The relations derived in reference 32 for predicting thermodynamic properties of 
real air in thermochemical equilibrium were also employed in the present study. These 
relations were obtained from curve fits to the real-air results of references 25, 33, 
and 34, and cover a temperature range of 90 K to 15 000 K. Imperfect air (intermolecular 
force) effects are neglected in reference 32. The maximum percentage errors in the 
results of these relations for a pressure range of 10 to 1000 kN/m2 and temperature range 


of 2000 K to 15 000 K are (ref. 32): 

a, percent 2.78 

h, percent 1.96 

T, percent 2.24 

Z, percent 0.75 

y E , percent =5.00 

p, percent 2.52 


These relations, in which the independent variables are p and s/R, are incorporated 
into a subroutine designated herein as SAVE (K). This subroutine utilizes an iteration- 
interpolation scheme allowing combinations of h, p, s/R, and p as inputs; K = 1 
denotes inputs p and s/R, K = 2 denotes inputs p and p, K = 3 denotes inputs p 
and h, and K = 4 denotes inputs h and p. 

Iterative Procedure for Standing Normal Shock 

The conservation relations for mass, momentum, and energy for a standing normal 
shock are 


PlVi = P 2 V 2 

(1) 

Pl+PlVi 2 = p 2 +P 2 v 2 2 

(2) 

h l4 V l 2==h 2 + i V 2 2 

(3) 
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Both the direct solution and the inverse solution to the conservation relations were used 
in this study. The direct solution requires the free-stream flow conditions, including 
flow velocity, as inputs and yields the postshock static and stagnation conditions. The 
inverse solution requires postshock flow conditions as inputs and yields the free-stream 
conditions. The direct solution is generally employed when at least two of the required 
three inputs are free-stream conditions, and the inverse solution is employed when at 
least two of the required three inputs are postshock conditions. Both solutions will now 
be discussed. 

Direct solution .- In the direct solution, the free-stream flow conditions appearing 
on the left side of equations ( 1 ) to ( 3 ) are considered known. To solve for the four 
unknown postshock static conditions (p 2 , V2, P2> and I12), an additional relation is 
required, the equation of state in the form of SEARCH (L) or SAVE (K). (Either 
SEARCH (L) or SAVE (K) may be used as the equation of state. In the subsequent dis- 
cussion SEARCH (L) will be used, except in cases where relatively accurate initial esti- 
mates of thermodynamic properties are required, for which SAVE (K) is preferable.) 

The iterative procedure for the direct solution begins by first estimating a value 
of P2. This value of P2 is used in equation ( 1 ) to yield a value of V2. Then pg is 
found from equation ( 2 ) and h2 from equation ( 3 ). These values of P2 and I12 are 
used as inputs to SEARCH ( 1 ) to obtain a value of P2. This P2 from SEARCH ( 1 ) is 
compared with the initial estimate of P2. If these values of P2 are not within 0.1 per- 
cent, the P2 obtained from SEARCH ( 1 ) is used in equations ( 1 ) to ( 3 ) to obtain new val- 
ues of V2, P2> and h2- This procedure, commonly referred to as the method of suc- 
cessive approximations, is repeated until successive values of P2 obtained from 
SEARCH ( 1 ) are within 0.1 percent. 

The number of iterations in the method of successive approximations is dependent 
upon the accuracy of the initial estimate of P2» Hence, a means for providing an accurate 
initial estimate of P2 was examined. Since the free-stream conditions are assumed to 
be known in the direct solution, the stagnation enthalpy can be determined from the 
relation 


h t = hi + |v 1 2 (4) 

If measured p^ is not available, the relation 

P t = CPjVi 2 (5) 
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is used to estimate p^. The factor C for real air may vary from approximately 0.93 to 
1.0 (ref. 35). In this study, C was taken to be 0.965. 

With a known p t value or estimated (within ±3.5 percent) value from equation (5) 
and a known fy, a relatively accurate value of p t is obtained from SAVE (3). The ini- 
tial estimate of P 2 is taken as 

P2 = 0.955P|. (6) 

In general, two to five iterations are required for this method of successive 
approximations. 

The postshock stagnation conditions are determined by assuming that the flow region 
extending from immediately behind the normal shock to the stagnation point is isentropic 
(that is, = S 2 /R). The calculated p^ is obtained by using s^R and ht as inputs 

to SLOW. Other stagnation-point thermodynamic properties of interest are obtained by 
using ht and p^. as inputs to SEARCH (1). 

For the case where SAVE (K) is used to obtain postshock stagnation conditions, a 
value of p^ is estimated from the ideal-gas isentropic relation (ref. 36): 

( y - 1 yEt/fax 1 ) 

Pt = p 2 ( 1+ - 1 T— m 2 2 | <?> 

Since the flow in the region extending from immediately downstream of the shock to the 
stagnation point can be considered to behave as an "ideal" gas (y^, 2 ~^E t anc * ^2 ~ ^t)> 

p^ from equation (7) is relatively accurate. Hence, values of p^. within 5 percent of p^ 
are used along with S 2 /R as inputs to SAVE (1), and corresponding values of h t are 
obtained. Since ht is known, an interpolation is performed to obtain the desired p^. 

For convenience, this iterative procedure for determining postshock static and 
stagnation-flow conditions will be referred to hereafter as DIRECT. 

Inverse solution .- The postshock conditions serving as inputs in the various proce- 
dures are inferred ht from stagnation- point heat-transfer rate measurement, measured 
Pt, and measured p^. For procedures having and p^. as inputs, the corresponding 
stagnation conditions are obtained from SEARCH (1). Procedures having p^ and as 
inputs utilize SEARCH (2) to obtain corresponding stagnation conditions. (Procedures 
having ht and p ^ as inputs are not included in this study.) 

The flow region behind the shock is assumed to be isentropic; hence S 2 /R is 
known. An initial estimate of P 2 is made and used with S 2 /R as input to SLOW to 
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obtain Pg and h 2 - The corresponding V 2 is found from the energy equation 


V 2 = f(ht ~ h 2 ) (8) 

In the data- reduction procedures using the inverse solution, either pj or Pj is known. 
Thus, the free-stream quantities appearing on the left side of equations (1) to (3) can be 
determined. The pj and hj values are used as inputs to SEARCH (1) to .obtain a value 
of Pj. This Pj from SEARCH (1) is compared with the measured or calculated pj, 
and if not within the desired tolerance, the value of P 2 is varied and the procedure 
repeated. This numerical iteration on P 2 is continued until the condition on Pj is 
satisfied. 

Because of the relative insensitivity of P2 to variations in free-stream flow condi- 
tions, this inverse shock- crossing procedure generally requires a larger number of itera- 
tions than the direct solution. Thus, it is desirable that the limits of P2 required in this 
iteration be chosen so as to minimize computer time. For most tests in the expansion 
tube, the limits on P2 will lie between 

(P 2 )up = °- 97 Pt 

(P2)low = 0-85Pt 


These limits are applicable for ideal-air free-stream Mach numbers greater than 3 or so. 
For Mach numbers less than 3, the lower limit of P 2 must be decreased to insure that 
the actual P 2 lies between these limits on pg. (In the interest of computer time, the 
user is urged to adjust these limits so as to minimize the range of iteration on P 2 for 
his particular problem.) 

This iteration procedure for determining postshock and free-stream flow conditions 
will be referred to hereafter as INVERSE. 

Prediction of Stagnation-Point Heat-Transfer Rates 

Over the past two decades, considerable effort has been directed toward obtaining 
relations for predicting stagnation-point heat-transfer rates at high flight velocities. 
Consequently, the literature contains numerous theoretical procedures for determining 
laminar stagnation-point heating rates for blunt axisymmetric bodies. However, as noted 
in references 37 and 38, the scatter in experimental data obtained at velocities greater 
than 6.1 km/sec prohibits identification of the least uncertain theoretical procedures. 
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Because measured q t is considered as a basic input datum, an accurate means of 
predicting must be utilized. Hence the theoretical results of references 39 to 43 for 

predicting on blunt axisymmetric bodies are examined, along with the empirical rela- 

tion of reference 44. (The results of these studies are discussed in appendix A.) As 
noted in reference 45, uncertainties in transport properties of high-temperature air 
represent a source of discrepancies in the various theoretical relations for predicting q t . 
Thus, the results of references 46 to 48 concerning transport properties of high- 
temperature air are discussed briefly in appendix B. 

Procedures for Determining Free-Stream and Postshock Flow Conditions 

The procedures for determining free-stream and post-normal-shock flow conditions 
are identified in the computer program as ITEST. For convenience, this method of iden- 
tification will be employed in the following discussion. 

The basic measured inputs and iterative procedure for crossing the normal shock 
are given in the following table: 


ITEST 

Measured inputs 

Shock- crossing procedure 

1 

P t ’ Pi’ V 1 

Direct 

2 

P t » V V 1 

Direct 

3 

Pt» %> p l 

Inverse 

4 

q t > v v p i 

Direct 

5 

p i’ p i’ v i 

Direct 

6 

p i> p i> p t 

Direct 

7 

p t> p t> p i 

Inverse 

8 

Pt» p t> p i 

Inverse 

9 

p t’ p t> v i 

Direct 


In these procedures calculated and measured flow quantities are compared, and if 
the calculated quantities are not within a prescribed tolerance of the measured quantities, 
a numerical iteration is performed. This iteration results in an upgrading of the cal- 
culated flow quantity until satisfactory agreement between calculated and measured val- 
ues is obtained. The prescribed tolerance established for iteration is, of course, some- 
what arbitrary. For example, in general usage of these procedures for the expansion 
tube, the calculated p t is required to be within 1 percent of the measured p^. Refining 
this tolerance to a smaller value is not believed warranted, considering the experimental 
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uncertainty in measured p^. However, in real-air facilities where pj. can be measured 
with a high degree of confidence, the user of these data- reduction procedures may wish to 
define the tolerances for iteration. Hence, the tolerances for iteration on the measured 
quantities p^., q^., and are treated as inputs and the user can establish his own 
degree of refinement on these iterations. The tolerance on measured p^ is denoted by 
TOLPT, that on measured q t by TOLQT, and that on measured p t by TOLRHO. 

In the case of the expansion tube, the free-stream static pressure is assumed to be 
equal to the measured tube wall pressure (that is, pj = p w This assumption is sub- 
ject to question. As discussed in reference 49, a number of experimental studies have 
shown that the measured wall static pressure is greater than the static pressure at the 
edge of a turbulent boundary layer. These studies indicate that p w ^jVi becomes 
increasingly greater than 1 as Mach number increases. For air and nitrogen flows, 

P w i/Pj is less than approximately 1.1 when is less than 10; however, when 

is about 20, p w ^ may be as much as twice p^ (ref. 49). Because of the inherent dif- 
ferences between conventional wind tunnels and the expansion tube (particularly in regard 
to test time), the questionable state of the expansion tube wall boundary layer (laminar, 
transitional, or turbulent), and the lack of a correlation for p w p^, no attempt is made 
to adjust p^. 

The calculation scheme for the individual procedures (ITEST) is discussed in 
appendix C. 


DISCUSSION 

Because of the greater simplicity associated with procedure ITEST = 5, this pro- 
cedure was the first to be programed for computer usage. For debugging purposes, cases 
were run with ITEST = 5 which covered a Tj. range of 1200 K to 14 000 K. The charts 
of reference 50 were used, where applicable, to provide a rough check on the computed 
flow parameters. After the debugging process for ITEST = 5, several cases were run 
and the results were compared with the recent results of reference 51. In all cases, the 
free-stream and postshock flow conditions agreed with those of reference 51 to within 
1 percent. (The Tt range corresponding to these comparison cases was 1500 K to 
13 000 K.) Following the successful check of ITEST = 5, the remaining eight procedures 
were programed and a common check case was run for each procedure. For these check 
cases the free-stream inputs were the same as those employed with ITEST = 5, but the 
postshock outputs p t , p^, and of ITEST = 5 were used as inputs. In the check 
cases of the various procedures, TOLPT = TOLQT = TOLRHO = 0.001. For these tol- 
erances of iteration, the computed flow quantities for all eight procedures were observed 
to be in excellent agreement with those of ITEST = 5. 
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The ITEST = 5 procedure was used to run a number of cases in which values of 
q^. were calculated from the theoretical findings of references 39 to 43 and the empirical 
result of reference 44. Included in these cases was a 15.24 km/sec entry trajectory for 
a vehicle having a lift-drag ratio of 1. The velocity range considered for this entry tra- 
jectory was 1.5 to 12.2 km/sec and the corresponding altitude range was 36.6 to 61 km. 

The results for this entry trajectory are shown in figure 2, where q t is plotted as a 
function of Vj. The value of T w was held at 300 K and r g (for a sphere) was 1.27 cm. 
The Tt value corresponding to Vj and to the thermodynamic conditions at the corre- 
sponding altitude can be found from the scale at the top of figure 2. Except for refer- 
ence 41, the source of transport properties used in the derivation of the various theoreti- 
cal results was also used for the predicted q t of figure 2. Instead of using Sutherland’s 
viscosity relation in the method of reference 41, was obtained from reference 46. 

At the lowest value of figure 2 there should be little, if any, difference in 
transport properties used in references 39 to 43. For this Vj the result of reference 44 
is approximately 1.2 times the average of the five theoretical results (q t ) av - As Vj 
increases, the difference between q^. from reference 44 and (q^ av decreases. The 
results of reference 46 are believed to be too low for T t values between 1500 K and 
8000 K (see appendix B); hence, if the more accurate ju t results of reference 47 are 
employed in references 39, 41, and 42, (q t ) av will increase. Between velocities of 4.5 
and 9 km/sec, the magnitude of this increase is such as to bring (q^ av into good agree- 
ment with the result of reference 44. At velocities above 9 km/sec (corresponding to 
Tt > 8000 K), the uncertainty in transport properties increases because of ionization phe- 
nomena. For values from 9 to 12 km/sec, the of reference 46 is still conserva- 
tive in comparison with that of reference 47; hence, use of the results of reference 47 
would increase ^t)av' T ^ s (^t)av wou ^d then be greater than the q^. of reference 44 
for this velocity range. 

In this study, the empirically based result of reference 44 for predicting q^. was 
decided upon for T^ = 4500 K. This corresponds to Vj greater than approximately 
4.5 km/sec in figure 2. (Although ref. 44 may be somewhat conservative for 
Vi > 9 km/sec, the uncertainty in transport properties at these conditions precludes 
modification of ref. 44 or adoption of another source for predicting q^.j For 
T t < 4500 K, the result of reference 41 ^with Nl 6 = 1^ was adopted, where is 
obtained from reference 47. 

The uncertainties in calculated flow quantities due to uncertainties in experimental 
inputs were examined. For each procedure, one of the inputs was varied (simulating a 
measurement error) while the remaining two inputs were held constant. Hence, the sen- 
sitivity of the calculated flow quantities to this variation in a given input was determined. 
The results of this error analysis are shown in figure 3, where the error parameter £ 
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is plotted for various free-stream quantities and stagnation-point quantities of interest. 
The parameter £ is defined as 


€ = 


error 

error 


^correct \ ^correct 
" ^correct/ ^correct 


0 ) 


where 6 represents the calculated flow quantity of interest and <p represents the 
experimental input quantity. The "correct" inputs used in obtaining the results of fig- 
ure 3 correspond to a representative expansion tube test with a heated helium driver at 
moderate driver pressure (16.5 MN/m2). (Values for the "correct" free-stream and 
postshock flow conditions are given in the sample printout of appendix D.) The terror 
was less than or equal to 5 percent for the results of figure 3. The tolerances of iteration 
used in obtaining these results were TOLPT = TOLRHO = 0.001 and TOLQT = 0.005. 

As observed from figure 3, the degree of uncertainty in calculated flow conditions 
corresponding to an uncertainty in experimental inputs varies for the different ITEST 
procedures. Naturally, the type of investigation being conducted in a facility would dictate 
what flow quantities are most important and what limits of uncertainty in these quantities 
can be tolerated. ^For example, Nr 6j i is an important parameter in most viscous flow 
studies but is of relatively little importance in stagnation-point radiative heating studies. j 
For purposes of illustration, let it be assumed that all the free-stream quantities shown 
on the abscissa of figure 3 are pertinent to a given investigation. The degrees of uncer- 
tainty in calculated free-stream quantities for the various procedures can be roughly 
compared by allowing the maximum uncertainty permissible in any of the calculated free- 
stream quantities of figure 3 to be some value, say 20 percent or so. For the represen- 
tative expansion tube test under consideration, the approximate maximum uncertainty 
permitted in each input (where the remaining two inputs are assumed to be correct) for 
each procedure is as follows: 


ITEST 

Maximum uncertainty, percent, in input 

— 

Pi 

V 1 

pt 

% 

Pi 

Pt 

1 

20 

8 

12 

— 

— 

— 

2 

— 

.5 

2 

1 

— 

— 

3 

20 

— 

7 

18 

— 

— 

4 

20 

2 

— 

5 

— 

— 

5 

20 

20 

— 

— 

12 

— 

6 

20 

— 

42 

— 

18 

— 

7 

— 

— 

2 

— 

1 

1 

8 

20 

— 

36 

— 

— 

10 

9 

— 

.5 

1 

— 

— 

1 
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Hence, for a study at the conditions given in the sample printout of appendix D, if all the 
free-stream quantities of figure 3 are considered pertinent, procedures ITEST = 1, 3, 5, 
6, and 8 are preferable to ITEST = 2, 4, 7, and 9. If the pertinent quantities for an 
investigation are the stagnation-point quantities of figure 3, the relative preferability of 
the procedures is not nearly so obvious. This is illustrated by the following table, where 
the maximum uncertainty permissible in any of the calculated stagnation-point quantities 
of figure 3 was taken as 20 percent or so: 


ITEST 

Maximum uncertainty, 

percent, 

in input — 


Pi 

V 1 

Pt 

<*t 

Pi 

Pt 

1 

»20 

10 

20 

— 

— 

— 

2 

— 

»5 

15 

35 

— 

— 

3 

»20 


15 

35 

— 

— 

4 

»20 

3.5 

— 

8 

— 

— 

5 

»20 

7 

— 

— 

20 

— 

6 

»20 


13 

— 

20 

— 

7 

— 


10 

— 

»20 

11 

8 

»20 


10 

— 

— 

12 

9 

— 

»5 

10 

— 

— 

11 


Error analyses were also performed for procedures ITEST = 1,3, and 5 for the 
same free-stream thermodynamic inputs that are shown in the sample printout of appen- 
dix D, but at stagnation conditions corresponding to V-^ values of 3 and 12 km/sec. 
(Procedures ITEST = 1, 3, and 5 were chosen because these are expected to be used 

0 error 

for most of the data reduction in the Langley 6-inch expansion tube.) The ratio 

^correct 

for both velocities (3 and 12 km/sec) was observed to be essentially the same as for the 
representative expansion tube test case considered previously for a velocity of 
6.1 km/sec. 

In a program such as that presented herein, computer time is of concern. The pro- 
cedures were run individually on a Control Data 6600 series computer for the represen- 
tative expansion tube test case used in the error analysis. The total time (computational 
and peripheral) for each procedure, when the thermodynamic properties were obtained 
from the AEDC tape (that is, SAVE (K) was not utilized), is given in the following table: 
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ITEST 

Total time, sec 

1 

580 

2 

490 

3 

375 

4 

570 (1200) 

5 

240 

6 

340 

7 

440 

8 

530 

9 

420 


These times were obtained with TOLPT = TOLQT = TOLRHO = 0.005 and the iteration 
limits presented herein. The total time of 570 seconds for ITEST = 4 corresponds to 
the refined limits on for greater than approximately 8; the total time of 

1200 seconds corresponds to the more general limits on Pj. (See appendix C.) 

Although these total times for procedures employing DIRECT should be fairly rep- 
resentative for the stated iteration tolerances, the time for procedures employing 
INVERSE (ITEST = 3, 7, and 8) should not. This is because P2/Pt = 0.959 for the exam- 
ple test case under consideration. Since the iteration procedure on P 2 in INVERSE 
begins with the upper limit on P 2 (^P2)up/Pt = 0-97^, the total times in the table above 

for ITEST = 3,7, and 8 are believed to be somewhat less than that corresponding to the 
general case. 

The relatively large computer times associated with the data- reduction procedures 
presented herein are due primarily to the time required for tape manipulation. This con- 
clusion was verified by examining the total time required by SEARCH (L). Cases were 
run for L = 1, 2, and 3. For a single usage of SEARCH (L), the total time required 
ranged from 34 seconds for L = 2 and 3 to 52 seconds for L = 1. These same cases 
were also run with SLOW. Multiple callings of SLOW were performed in order to obtain 
the same thermodynamic properties as were obtained with SEARCH (L). The total time 
was 28 seconds. Hence, it is obvious that repetitive usage of SEARCH and SLOW, as 
required in the iterations of the present procedures, will consume a large amount of total 
computer time. 

Subroutine SAVE (K), which is based on the real-air curve-fit expressions of refer- 
ence 32, was incorporated into the present data-reduction procedures in an attempt to 
reduce computer time and tape usage. This subroutine replaces the subroutines SLOW 
and SEARCH (L), which were written to search the AEDC real-air tape, as the source 
of real-air thermodynamic properties. Some loss in accuracy is incurred in using 
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SAVE (K), particularly at high densities where inter molecular force effects become 
important. (See section entitled "Thermodynamic Properties for Real Air.") However, 
these data- reduction procedures based on SAVE (K) should be sufficiently accurate for 
most purposes. The nine procedures, now divorced from the AEDC tape, were run indi- 
vidually on the computer for the same representative expansion tube case and same iter- 
ation tolerances used in the previous time study. For each procedure the computer time 
was much less than when the AEDC tape was used. For example, the 1200 seconds 
required for ITEST = 4 with subroutines SLOW and SEARCH (L) was reduced to 40 sec- 
onds with subroutine SAVE (K). Thus, the problem of relatively large computer times is 
circumvented by usage of SAVE (K) without sacrificing appreciable accuracy. 

CONCLUDING REMARKS 

Data- reduction procedures for determining free-stream and post-normal-shock 
kinetic and thermodynamic quantities are derived. These procedures are applicable to 
imperfect real-air flows in thermochemical equilibrium for temperatures to 15 000 K and 
a range of pressures from 0.25 N/m2 to 1 GN/m2. Nine data-reduction procedures 
resulting from various combinations of three of the following measured flow parameters 
were derived: 

(1) Stagnation pressure behind normal shock 

(2) Free-stream static pressure 

(3) Stagnation-point heat-transfer rate 

(4) Free-stream flow velocity 

(5) Stagnation density behind normal shock 

(6) Free-stream density 

The various combinations of measured flow parameters are identified herein as ITEST 
and are: 


ITEST 

Measured flow parameter 

1 

(1), (2), (4) 

2 

(1), (3), (4) 

3 

(1), (2), (3) 

4 

(2), (3), (4) 

5 

(2), (4), (6) 

6 

(1), (2), (6) 

7 

(1), (5), (6) 

8 

(1), (2), (5) 

9 

(1), (4), (5) 
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These procedures employ an adjustment of computed flow parameters by numerical 
iteration until measured and computed flow parameters are within a prescribed toler- 
ance. All nine procedures are incorporated into a single computer program written in 
FORTRAN IV language. 

The uncertainties in calculated flow quantities due to uncertainties in the experi- 
mental inputs were examined. This error analysis demonstrated that for an investiga- 
tion in which free -stream quantities (including Reynolds number) were pertinent, pro- 
cedures ITEST = 1, 3, 5, 6, and 8 are preferable to ITEST = 2, 4, 7, and 9. 

Relatively large computer times were observed for these procedures. The large 
times are due, primarily, to the time required in searching the real-air tape. Significant 
reduction in computer time, without appreciable loss of accuracy, was obtained by using 
real-air curve-fit expressions as the source of thermodynamic properties in place of the 
tape. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 11, 1972. 
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APPENDIX A 


RELATIONS FOR PREDICTING STAGNATION- POINT HEAT-TRANSFER RATES 


In the present study the theoretically derived expressions of references 39 to 43 for 
predicting stagnation-point heat-transfer rates on blunt axisymmetric bodies were exam- 
ined. The results of these five theoretical studies are discussed in this appendix, and the 
empirical stagnation-point heat-transfer relation of reference 44 is presented. 


Cohen 

From correlations of numerical results, Cohen (ref. 39) obtained the relation 



for predicting stagnation-point heat-transfer rate. Cohen considered two free-stream 
velocity regimes in the derivation of equation (10). In the lower velocity regime, where 
the velocity was less than 8.84 km/sec, the air was assumed to be an equilibrium mixture 
of oxygen and nitrogen atoms and molecules. The transport properties for this equilib- 
rium dissociated air were taken from the correlations of reference 52. In this lower 
velocity regime the results were obtained for wall temperatures from 300 K to 1750 K, 
and equation (10) represents the numerical solutions to within approximately ±5 percent. 
In the upper velocity range, where the velocity was from 8.84 to 12.5 km/sec, the trans- 
port properties of Hansen (ref. 46) were used by Cohen. Results were obtained for wall 
temperatures to 5200 K, and equation (10) represents the numerical solutions for this 
regime to within approximately ±10 percent. 

In the present study, the thermodynamic quantities and appearing in equa- 
tion (10) are obtained from the AEDC real-air tape. The value of ju.^ is obtained from 
reference 46 or 47. (Further discussion concerning the obtainment of is presented 
in appendix B.) The parameter e appearing in equation (10) distinguishes the velocity 
regime. For the lower velocity regime (to 8.84 km/sec), e = 0; for the upper velocity 
regime (8.84 to 12.5 km/sec), e = 1. 

Because of the very short test time of the expansion tube, the wall temperature 
remains on the order of ambient temperature. Hence, the wall temperature can be set 
equal to the ambient temperature, as was done in reference 22, or can be estimated. To 
obtain a rough estimate of T w , the assumption is made that the thick- film heat-transfer 
gage used to measure stagnation-point heat-transfer rate experiences uniform flow. This 
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APPENDIX A — Continued 

results in a constant heat flux during the test time. Then, from reference 17, the wall 
temperature is given by the expression 


T w 



+ lamb 


( 11 ) 


(This expression was derived from the one -dimensional heat- conduction equation for a 
homogeneous semi -infinite slab.) 

Since T w remains on the order of T am ^, the viscosity at the wall is obtained 
from Sutherland's expression (ref. 46): 




1.462 x 10" 6 



(12) 


The density at the wall is determined from the equation of state, 


P t W 0 

P Wrt^ 

and the static wall enthalpy from the ideal-air expression, 


(13) 


h w - 3.5 T w (14) 

It should be noted that in most cases hf. » h w , and too, the dependence of equation (10) on 
p w jn w is to the 0.07 power. Hence, the stagnation- point heat-transfer rate obtained from 
equation (10) has a weak dependency on wall temperature. For the range of Tw expected 
in expansion tube or expansion tunnel testing, relatively large errors in Tw would not be 
expected to have a significant effect on the calculated stagnation-point heat-transfer rate. 

The stagnation-point velocity gradient was calculated from the modified Newtonian 
relation 



(15) 


for axisymmetric bodies having a nose radius r^. 

Because of the relatively small range of Tw expected in the expansion tube or tun- 
nel, the Prandtl number at the wall was assumed to be a constant and was set equal to 0.71. 
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APPENDIX A — Continued 


Hoshizaki 

Hoshizaki (ref. 40) correlated his numerical results to within ±6 percent, deriving 
the relation 


q t = 1.1672\|^p w p. w T w °’ 4 V 1 1 ’ 69 |l - (16) 

According to reference 40, this relation is valid for velocities from 1.83 to 15.2 km/sec, 
stagnation pressures between 100 N/m2 and 10 MN/m^, and wall temperatures from 300 K 
to 3000 K. The boundary layer at the stagnation point is assumed to be in thermochemical 
equilibrium. The effects of dissociation and ionization are taken into account by use of the 
total thermodynamic and transport property concept. In this concept, the properties are 
defined in such a way that the effects of dissociation and ionization are contained within 
them. The calculations of Hoshizaki are based on the transport properties of Hansen 
(ref. 46) for viscosity and allow the Prandtl number and Lewis number to vary. Methods 
employed herein for determining the various quantities appearing in equation (16) were 
discussed in the preceding section. 


Fay and Riddell 


Fay and Riddell (ref. 41) obtained, from correlation of numerical results, the 
expression 


q^ - 0.76^Np r>w ) (Pt^t) (P w Pw) (hfh w )^ 


1 + N 


J Le 


0.52 



(17) 


for predicting stagnation-point heat-transfer rate for dissociated air in thermochemical 
equilibrium. These results were restricted to a wall Prandtl number of 0.71; wall Lewis 
numbers of 1.0, 1.4, and 2.0; wall temperatures from 300 K to 3000 K; and velocities 
from 1.77 to 6.95 km/sec. The viscosity used in reference 41 was obtained by using 
Sutherland's expression (eq. (12)). 

In reference 38, a semiempirical modification of Fay and Riddell's theory is sug- 
gested to extend it to higher velocities where the effects of ionization become important. 
Lewis and Burgess (ref. 38) compared the theory of Fay and Riddell for equilibrium nitro- 
gen (for nitrogen, the factor 0.76 in eq. (17) is replaced by 0.754) and N Le = 1.0 with the 
results of Fay and Kemp (ref. 53). Fay and Kemp utilized a simplified binary diffusion 
model of an ionized diatomic gas to obtain the transport properties for nitrogen. The cal- 
culations of reference 53 were performed for a wall temperature of 300 K, and a Lewis 
number of 0.6 was employed in most of the calculations. The comparison of Lewis and 
Burgess (ref. 38) showed that if the results of Fay and Riddell are simply multiplied 
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APPENDIX A - Continued 


by 1.15, they are in good agreement with the theory of Fay and Kemp (which includes 
ionization effects) for velocities in nitrogen from 6.1 to 16.78 km/sec. However, on the 
basis of the summary plot of reference 38, the present authors do not believe such a 
semiempirical modification to the theory of Fay and Riddell is justified for this study. 
Instead, it is felt that the theory of Fay and Riddell for equilibrium air and a Lewis num- 
ber of 1 provides a "compromise" means for estimating the stagnation-point heat-transfer 
rate for velocities to 13.7 km/sec. (This conclusion is based on a comparison of curve 1 
with the experimental data in the summary figure of ref. 38.) 

The quantities appearing in equation (17) for Nl 6 = 1 are determined in the same 
manner as those in equation (10) derived by Cohen. 

Pallone and Van Tassell 

From a correlation of numerical results, Pallone and Van Tassell (ref. 42) derived 
the relation 


q t = 0.90(Np r , W )’ 0 ' 75 (P t M t )°‘ 43 (P w M w ^ 0 ‘ 07 (ht - 
where for V^ = 9.906 km/sec, A = 1, and for Vi > 9.906 

a.AL 

9906 

The parameter A represents a velocity dependence of the boundary -layer solutions for 
velocities in excess of 9.906 km/sec. Above velocities of 12.5 km/sec, a slight depen- 
dence on stagnation- point pressure was also observed. The calculations of Pallone and 
Van Tassell for air are also based on the transport properties of Hansen (ref. 46). 
Methods employed for determining the various quantities appearing in equation (18) were 
discussed previously. 


h w )f§A (18) 

km/sec, 


DeRienzo and Pallone 

The primary objective of this more recent study by DeRienzo and Pallone (ref. 43) 
was to extend the results of reference 42 to flight speeds as high as 21.3 km/sec. Unlike 
the studies of references 39, 40, and 42, this study (ref. 43) utilized the more recent 
transport properties of reference 47. Calculations were performed for velocities from 
1.52 to 21.3 km/sec and included the effects of blowing. The numerical results of refer- 
ence 43 are given in table I of that reference. To obtain expressions for calculating the 
stagnation-point heating rate, these tabulated results were correlated by the present 
authors for the case of an axisymmetric body with no blowing, stagnation-point pressures 
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of 0.1 and 1 MN/m^, and a velocity range of 1.52 to 12.2 km/sec. The results of refer- 
ence 43 were correlated to within 5 percent by the expression 


n Nu,w 

^ N Re,w,S 


0.62 


'Pt*V 

,P w Mw, 


0.3376 


for the heat-transfer parameter. The basic relation for qj- is (ref. 54). 


(19) 


. T -l n Nu,w i -/, , \ 


( 20 ) 


and substituting equation (19) into equation (20) yields the expression 


q, = 0.621^^ .) 0 - 3376 (P w (x w )°' 1624 (h, - K)f (21) 

For velocities from 12.2 to 21.3 km/sec, the results of reference 43 were correlated to 
within 5 percent and the heat-transfer parameter was found to be 

%u,w 
\/ N Re,w,S 

The value of in equations (21) and (22) is obtained from a numerical interpolation of 
the viscosity data of reference 47. 


¥t 




(22) 


Zoby 

In reference 44, Zoby presents a simple empirical relation for predicting 
stagnation- point heat-transfer rates in several gas mixtures, including air. This rela- 
tion has the form 



where r e is the effective nose radius and Kj is a constant intended to account for the 
effect of thermodynamic and transport properties of the gas at the wall and external to 
the boundary layer. The constant Ki was determined in reference 44 by fairing a 
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straight line through results of prediction methods and experimental measurements on a 
plot of q^r e y/p t as a function of h^. - h w . 

The effective nose radius is defined in reference 55 as the hemispherical radius 
which produces the same velocity gradient as that computed for a blunt body which is not 
hemispherical. Since hemispherical heat-transfer probes will be used in the Langley 
6-inch expansion tube and expansion tunnel, the effective radius is identical to the geo- 
metric nose radius ^that is, r e = r^j. 

Transverse cylindrical heat-transfer probes will also be used. For these probes 
(assuming the flow about the cylinder is two dimensional), the stagnation-point heat- 
transfer rate is related to that of a sphere by the expression (from ref. 56) 

(^t) sphere " /^(^t)cylinder 

For air, equation (23) takes the form 



APPENDIX B 


TRANSPORT PROPERTIES FOR HIGH-TEMPERATURE AIR 

As discussed in reference 45, discrepancies in the results of theoretical stagnation- 
point heat-transfer studies can generally be attributed to (1) the assumptions employed to 
reduce the governing equations to a tractable form, (2) the mathematical technique used 
to solve the governing equations, and (3) the source of thermodynamic and transport prop- 
erties used. Aspects (1) and (2) have been fairly well standardized. However, as pointed 
out in reference 45, uncertainties exist in the evaluation of transport properties of high- 
temperature (ionized) air and may result in appreciable uncertainty (up to 20 percent or 
so) in the calculated stagnation-point heating rate. (It should be noted that the findings 
of ref. 57 indicate a more insensitive relation between transport properties at high tem- 
peratures and calculated q^ than the findings of ref. 45.) 

Cohen (ref. 39) states that as more accurate values of transport properties for real 
air become available, a better estimate of stagnation-point heat-transfer rate should be 
possible with equation (10). Cohen suggests that the coefficient and exponents appearing 
in his expression (eq. (10)) should be valid with these more accurate transport properties. 

As shown in the recent study of reference 48, the viscosity data of Hansen for tem- 
peratures above 1500 K appear to be too low, deviating from the results of reference 48 
by as much as 25 percent for temperatures to 8000 K. The viscosity results of refer- 
ence 47 agree closely with those of reference 48. In reference 48, where dissociation but 
not ionization phenomena are considered, this discrepancy is attributed to the fact that 
Hansen uses simple kinetic theory and rough approximations for collision cross sections. 

Because of the inexact knowledge of the transport properties of high-temperature 
air and the belief of Cohen that existing expressions may be valid with usage of more 
accurate transport properties, the viscosity results of both reference 46 and reference 47 
are made available to the expressions of references 39 to 43. This permits the user of 
the procedures of this study to choose transport properties from reference 46 or 47 in 
predicting q^. (see appendix E). 
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APPENDIX C 


CALCULATION SCHEME FOR INDIVIDUAL DATA-REDUC TION PROCEDURES 

ITEST = 1 

In the ITEST = 1 procedure, pj and Vj are measured inputs, and a second 
free-stream thermodynamic quantity is determined. This is accomplished by estimating 
a value of Pj from equation (5), since p^. is also a measured input. This value of Pi 
is believed to be within approximately ±3.5 percent of the actual value. Then p^ and pq 
are used as inputs to SAVE (2) to obtain the corresponding free-stream thermodynamic 
quantities. 

With known values of pj and Vj and relatively accurate values of Pj and hj, 
the DIRECT iterative procedure is performed. The calculated p^. obtained from 
DIRECT is compared with the measured p^. and if not within TOLPT, h^ is adjusted 
by using the relation 


_ (^l)prev(^t)c,prev 
1 ’ ( p t)m 


( 26 ) 


This new value of hj is used, in conjunction with pj_, as input in SEARCH (1). The 
corresponding value of p^ obtained from SEARCH (1) is used in DIRECT. ^Note that 
p^ from eq. (5) was used only to obtain a relatively accurate first estimate of hj and 
is not involved in the final phase of upgrading hj.^ The from DIRECT is again 

compared with (p^ m - This procedure of varying hj according to equation (26) is con- 
tinued until is within TOLPT of ^P^ m * 

Additional free-stream parameters of interest are aj, s-j^R, Tj_, Zj, Z^*, and 
^E i’ These parameters are all included on the AEDC tape and thus are available from 

SEARCH (1) for the final value of hj and the known p^. Other free-stream parameters 
of interest are Mach number, coefficient of viscosity, and Reynolds number. The Mach 
number is found by dividing V^ by a^. For Tj = 1500 K, the coefficient of viscosity 
p is obtained from Sutherland's expression (eq. (12)), whereas for Tj > 1500 K, p is 
obtained from interpolation of the results of reference 47. The unit Reynolds number is 
found by dividing the product PjVj^ by pj. 
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ITEST = 2 


In the ITEST = 2 procedure, the only known free-stream flow quantity is Vj and 
the known postshock conditions are p^ and q^.. 

Since the procedure ITEST = 5 (to be discussed subsequently) was the first to be 
programed for computer usage, it was used to calculate q t from the theoretical expres- 
sions of references 39 to 43 and the empirical expression of reference 44. (See appen- 
dix A.) The five theoretical predictions, the average of these predictions, and the empir- 
ical prediction of Zoby were compared. (See fig. 2.) This comparison led to the adoption 
of the empirical relation of Zoby (eq. (25)) for Tt = 4500 K. For Tt < 4500 K, the the- 
oretical expression of reference 41 (eq. (17)) was adopted, with Nj_, e = 1 and with pt 
obtained from Sutherland's viscosity expression (eq. (12)) for Tj- = 1500 K and from 
interpolation of the results of reference 47 for 1500 K < T^. < 4500 K. 

The calculation scheme for ITEST = 2 begins by determining ht from equa- 
tion (25), with p^, q t , and r^ as measured inputs and h w obtained as discussed in 
appendix A. Then and pt are used as inputs to SEARCH (1) to obtain the postshock 
stagnation conditions. 


If Tf. = 4500 K, these stagnation conditions from SEARCH (1) are assumed to be the 
correct values. If T^ < 4500 K, the relation of reference 41 (eq. (17)) is used to obtain a 
value of h^, as will be discussed subsequently. The postshock stagnation conditions that 
have been found and the known value of Vj are used to obtain a value of hj from equa- 
tion (4). An estimate of pj is obtained from equation (5) and used in conjunction with 
hj as input to SEARCH (3) in order to obtain the corresponding free-stream conditions. 
The DIRECT iterative procedure is used to determine the postshock static conditions and 
a value of (Pt) c * If this (Pj) c is not within the prescribed tolerance (TOLPT) of the 
measured p^., the free-stream static pressure is adjusted according to 


( p l)prev( p t)m 

Pl=— 7-T (27) 

( p tjc,prev 

This value of p-^ and the previously determined hj are used with SEARCH (1) to obtain 
upgraded values of free-stream flow quantities, which are then used in DIRECT. This 
iterative procedure, based on equation (27), is continued until the desired agreement 
between (p^ m and (P^ c is obtained. Additional free-stream and postshock conditions 

of interest are determined as for ITEST = 1. 

For the case where T t < 4500 K, initial estimates of lit and Pt are obtained by 
using the ht value from equation (25) and the measured pt as inputs to SEARCH (1) to 
obtain Pf Then p t is obtained from equation (12) for Tt = 1500 K and from the 
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results of reference 47 for 1500 K < T^- < 4500 K. The /lx w , p w , and h w values are 
found from equations (12), (13), and (14), respectively. The value of is initially set 
at zero in equation (15). Since in most instances p^. » p-^ (for example, p^. > 10p^ 
for > 2.75 in ideal air), neglecting Pj in equation (15) will not result in an appreci- 
able error. (The term p t - pj appears to the 0.25 power in eq. (17).) The calculated 
q^. from equation (17) is compared with the measured q^. If not within the desired tol- 
erance, hj- is varied according to the relation 


ht = 


(ht - hw)prev(^t)m 
(«t)c 


+ h w 


(28) 


and the procedure is repeated. After is obtained, the free-stream conditions are 
determined as discussed previously. A value of /3 is calculated with pj included and 
compared with the value of /3 when p^ = 0. If these /3 values are not within a pre- 
scribed tolerance, the upgraded value of /3 is used in equation (17) and the procedure is 
repeated. 


ITEST = 3 

The measured postshock conditions for the ITEST = 3 procedure are p^. and q^. 
Hence the method for obtaining the postshock stagnation conditions is the same as that for 
ITEST = 2. The INVERSE iterative procedure is then used. In INVERSE, a value of P 2 
is estimated and the corresponding postshock static conditions are obtained. The mea- 
sured pj is used to obtain Vj by combining equations (1) and (2) to yield 


Vl 


P2 “Pi 

P 2 V 2 


+ V 2 


(29) 


Then p^ is found from equation (1) and h^ from equation (3). This h^ and the mea- 
sured p^ are used as inputs to SEARCH (1) to obtain the corresponding free-stream 
conditions. If the p-^ value from SEARCH (1) is not within the desired tolerance 
(TOLRHO) of Pi obtained from equation (1), p£ is varied and the procedure repeated. 


The iterative procedure on pg in INVERSE is as follows: An upper and a lower 
limit on P 2 are established as discussed previously. The upper limit is taken as the 


Ap“ = 


P2 

(that is, 

P 2 “ = 



f n \Ot 

a _ 

( p 2)up “ 1 

^2/low 


(30) 
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With pg**, a value of Pj is obtained as discussed above. If 


( p l)eq. (1) 
( P l) SEARCH (1) 


S TOLRHO 


( 31 ) 


then the value of p 2 a is considered satisfactory. If the condition of equation (31) is not 
met, and if SEARCH (1) > ( P l)eq. (1)’ then 


„ ct + 1 „ a A „ar 

P 2 = P 2 " A P 


(32) 


but if (Pi)seARCH (1) < ( P l)eq (1)’ the limits of p 2 are varied according to 




a+1 a 
ow " p 2 


(33a) 


(P2)up +1 = p 2 a + A P a! 


(33b) 


A new value of Ap is calculated ^Ap a+ ^), and P2^ + * becomes 

P 2 “ +1 - (P 2 )np +1 - ^“ +1 ( 3 “> 

ru i 1 

The procedure is repeated with p 2 and new values of are obtained. This iter- 
ative procedure is continued until the p^ value obtained from equation (1) is within the 
desired accuracy of the Pj value from SEARCH (1). 

In the course of this iteration, it is possible that values of hj less than the mini- 
mum value of hj on the AEDC tape (or even negative values of h^ may occur. In this 
case, the p 2 for which this occurred is varied according to equation (34) and the calcu- 
lation scheme is continued. 


ITEST = 4 

The measured inputs for the ITEST = 4 procedure are p^, Vi, and q t . The 
first consideration is to obtain an estimate of hi or p-^. Combining equations (4) 
and (25) and solving for p^. gives 
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r g (2. 57745 X 10 3 q t ) 2 


( 35 ) 


Assuming 


I>1 - « \ Vi 2 


and introducing equation (5) into equation (35) yields 


2.6841 x 10 7 r_q, 2 

Pj= — (36) 

Vl 6 

^An expression for hj in terms of the measured inputs was obtained by combining equa- 
tions (4), (5), and (25). However, this expression proved to be unsatisfactory since rela- 
tively small errors in q^. and Vj often resulted in hj < O.j The uncertainty in p^ 

1 9 

is dictated by the validity of the assumption hj - h w « — Vi , by whether Tj. = 4500 K 

z 

(region where eq. (25) is considered valid), and by the relatively small uncertainty in C 

1 9 

of equation (5). Calculations for typical expansion tube tests, where hj - h w « — V\* 

z 

and Tf- = 4500 K showed that Pj was approximately 2 to 5 percent greater than the 

actual value. However, calculations for a few shock tube cases, where T(- § 4500 K but 

hj - h w was not much less than — V^, showed that p^ was much smaller than the 

2 

actual value. For this procedure, limits on p^ that should prove to be valid for nearly 
all cases are 


( p l)up = 2 - 05? l 
( p l)low = °- 20p l 

For Mach numbers greater than 8 or so, these limits may be refined to 
( p l)up = 1>05p l 
( p l)low = °- 85p i 
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Again, the user is urged to adjust these limits so as to minimize the range of iteration 
on Pj for his particular problem. 

The Pj and measured p^ are used as inputs to SEARCH (2) to obtain the corre- 
sponding free-stream conditions. With estimates of Pj and hj and known values of 
Pj and Vj, the DIRECT iterative procedure is used to find postshock flow conditions, 
including q^.. This calculated q t is compared with the measured q t , and if not within 
the desired tolerance (TOLQT), the value of Pj is varied and the procedure repeated. 

The iterative procedure on Pj for the case where >8 is as follows: When 
the upper and lower limits on Pj have been established, the initial value of Pj is taken 

to be 0.95pj (that is, p^ a = 0.95p\ A value of A p a is obtained from 


Ap a = 



(37) 


OL 

With Pj and p^, the corresponding free-stream conditions are obtained from 
SEARCH (2), and DIRECT is used to find the postshock conditions. A value of q^ is 
calculated. If 


1 - 



^ TOLQT 


(38) 


then the value of p^ a is considered satisfactory. If the condition of equation (38) is not 
met, and if (q t ) c > (q t ) m , then 

p,“ +l = p,“ - 4p“ (39) 


but if (q t ) c < (q t ) m , the limits of 


pj become 



a+1 

up 



+ A p 


a 


(40a) 



Q!+l 

low 


= P 


a 

1 


A new value of A p is calculated 



and 



becomes 


(40b) 
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Pl“ +1 - (Pl)up +1 - ^“ +1 («« 

This iterative procedure is continued until the condition of equation (38) is satisfied. 

ITEST = 5 

The inputs for the ITEST = 5 procedure are p^, p^, and Vj. The values of p^ 
and Pj are used as inputs to SEARCH (2) to obtain the corresponding free-stream con- . 
ditions. Then the postshock conditions are obtained from DIRECT. 

ITEST = 6 

In the ITEST = 6 procedure, the measured inputs are p^, p^, and p^.. The 
free-stream thermodynamic conditions are obtained as for ITEST = 5. The value of Vj_ 
is estimated from equation (5) and the DIRECT iterative procedure is employed. The 
calculated p^. from DIRECT is compared with the measured p^. and, if within the 
desired tolerance (TOLPT), the estimated Vj is considered satisfactory. If not within 
the desired tolerance, the Vj value is adjusted by using the formula 



(42) 


This new value of is used in DIRECT. Then Vj_ is upgraded according to equa- 
tion (42) until is within TOLPT of (p^ m » 

ITEST = 7 

The measured inputs for the ITEST = 7 procedure are p-^, p^., and p^. The 
postshock stagnation conditions are obtained by using p^. and Pj. as inputs to 
SEARCH (2). The INVERSE iterative procedure is employed. When P 2 has been esti- 
mated and the corresponding postshock static conditions have been determined, is 
obtained from equation (1), pj from equation (2), and hj from equation (3). Then p^ 
and hj are used as inputs to SEARCH (1), and the corresponding p^ from SEARCH (1) 
is compared with the measured Pj. If not within the desired tolerance (TOLRHO), pg 
is varied and the procedure repeated. The variation of pg is the same as in the 
ITEST = 3 procedure, with the condition (P]J C < (Pj) m replacing the condition 

(^SEARCH (1) < (^eq (1) ITE ^T =3. In the course of this iteration on P 2 , it is 
possible that negative values of p^ and hj may occur. If so, then the limits on P 2 
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are adjusted according to equation (32), a new value of Pg is found, and the calculation 
scheme is continued. 

ITEST = 8 

The measured inputs for the ITEST = 8 procedure are p^, p^, and p t . The 
postshock stagnation conditions are determined as for ITEST = 7. The INVERSE itera- 
tive procedure is employed in the same way as in ITEST = 3. 

ITEST = 9 

In the ITEST = 9 procedure, the measured inputs are Vj, p t , and p t . The 
postshock stagnation conditions are determined as for ITEST = 7. Then hi is obtained 
from equation (4) and Pi is estimated from equation (5). The corresponding free- 
stream conditions are obtained from SEARCH (3). The DIRECT iterative procedure is 
employed and pj is upgraded as in ITEST = 2 (that is, according to eq. (27)) until 
(p t ) c is within the desired tolerance of 
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LISTING OF COMPUTER PROGRAM FOR DATA-REDUC TION PROCEDURES 

WITH SAMPLE DATA PRINTOUT 

The data-reduction procedures for determining free-stream and post-normal-shock 
flow conditions for real air in thermochemical equilibrium are incorporated into a single 
computer program. This program is written in FORTRAN IV language for a Control Data 
6600 series computer. Minimum machine requirements are 130 000 octal locations of 
core storage. A listing of this program, including subroutines and comments, is repro- 
duced on the following pages. 
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CO 

co 


JOB . 1 <1000, 130000 <6000* A3238 RGK 1 43 1247A CENT 

USER. MILLER. CHARLES G ITT 0006Q5575N 3454Q 

Rl IN ( S ) , 

RE QUEST tTADE8iHY« 70sQ42 . ROL . CGM . 

RFWIND(TAPfB ) » 

SET T NDF « . 

LGO. < 

UNLOAD ( T APP8 ) . ;J r 

PXlT.t 

UNLOAD ( TAPP8 ) . " 


PROGRAM M I LLER ( I NPUT . OUTPUT < TAPE5= I NPUT . T A PE 6 = OUTPUT . TAPE 8 ) A 1 

COMMON /BLK1 / RHQ5 . He , S5R < T5 . A5 . Z5 . GAME5 . ZST AR5 , I CODE A 2 

COMMON /BLK2/ RH0T5.HT5.ST5R.TT5.AT5.ZT5.GAMET5.ZSTART5.PT5 A 3 

COMMON /BLK3/ TWE . TAm TW < QT5DE . SPH < TWALL . ETA . SAV ; A 4 

COMMON /BLK4/ NV < IT « R < MM . I SP . NC O . MN . GAM E5S A 5 

COMMON /BLK5/ P5.U5.M5.MU5.RE5.QT5FR.QT5H0.S.QT5C0.QT5Z0. QT5PT A 6 

COMMON- /BLK6/ P5S . T5S . RH05S . H5S < A55 . Z5S . U5S . M5S . MU5S . RE5S . S5S R A 7 

COMMON /BLK7/ RN iTOLPT tTOLQT . TOLRHO tYOSt QD O . I TEST A 8 

DIMENSION RESULT (2) ‘ A 9 

REAL M5<MU5<M5S<MU5S . MUW i MUT5 ♦ NU .LAM A 1 0 

NAMELIST /INP/ P5M , USM i RH05M . P T 5M , QT5M . RH0T5M . I TEST . RN . TWE < T AU . TOL A 11 

1 PTtT OLQT <TOLRHO<YOS 1 ODO< SPH 1 TWALL 1 ETA 1 RIJN 1 .SAV A I 2 

CALL DAYTIM (RESULT) A 1 3 

1 P5M=PT5M=RH05M=RHOT5M=QT5M=U5M=O.Q A 1 4 


tolpt=tolrho=.oqi 
TOLQT=< 005 

sph=twall=o . 

YQS=ODO=SAV= 1 . 

TWF = gQO « 

RN= » O 1 

FTA = 3.043'' + 8 


TAU=.00Q1 



REA D ( 5 ♦ INP ) 


A 

15 


IF (FNDFILF 5) 16,2 


A 

16 

p 

ICODP=0 


A 

17 


NV “9 


A 

1 8 


I T = 8 


A 

19 


R = 2B'7*02^F 


A 

20 


MM =0 


A 

21 


ISP=1 


A 

22 


NCO = n 


A 

23 


PRINT 17. RFSUL T ( 1 ) 


A 

24 


PRINT 18 


A 

25 
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PRINT 19 



A 

26 


PRINT 20 



A 

27 


PRINT 21 



A 

28 


PRINT 22 ♦ RUN * P5M ,U5M « PT5M , QT5M . RH05M , RH0T5M , TAU 


A 

29 


GO TO (3«4*4*7*8*8« 1 1 « 1 1 * 1 1 > * ITEST 



A 

30 

r 




A 

31 

r 

ITEST 1 contained IN MAIN program 



A 

32 

c 




A 

33 

3 

PT5=PT5M 



A 

34 


P5 =P^M 



A 

35 


U5 = U c; M 



A 

36 


R5EST=PTF/( .965*U5**?) 



A 

37 


RH05=R5EST 



A 

38 

c 




A 

39 

C 

IF SAv=0* THERMODYNAMIC PROPERTIES OBTAINED 

FROM 

AEDC TAPE 

A 

40 

c 

IF S A V= 1 * THERMODYNAMIC PROPERTIES COMPUTED 

FROM 

CURVE-FIT 

A 

41 

c 




A 

42 


CALL SAVE (P5.RH05.H=s,S5R,T5.A5»Z5iGAME5i2) 



A 

43 


CALL DIRECT (PT5M ) 



A 

44 


GO TO 15 



A 

45 

4 

IF flTEST.EO.3) GO To 5 



A 

46 


XX=U5M 



A 

47 


GO TO 6 



A 

48 

c; 

XX = P«=:M 



A 

49 

6 

CALL PR0C2 (PT5M,QT5m,XX) 



A 

50 


GO TO 1 5 



A 

51 

7 

CALL PROC4 (P5M »U^M ,OT5M ) 



A 

52 


GO TO 15 



A 

53 

8 

IF (ITEST. EQ. 5) GO To 9 



A 

54 


YY = PT5M 



A 

55 


GO TO 10 



A 

56 

Q 

yy=u*m 



A 

57 

10 

CALL PROC5 (RH05M*P5m«YY) 



A 

58 


GO TO 15 



A 

59 

•1 1 

IF (ITEST. EQ.7) GO TO 12 



A 

60 


IF (ITEST. EO. 8) GO To 13 



A 

61 


YY = IJ5M 



A 

62 


GO TO 14 



A 

63 

1 2 

YY=RH05M 



A 

64 


GO TO 14 



A 

65 

1 3 

YY=P5M 



A 

66 

14 

CALL PR0C7 (PT5M,RH0T5M,YY ) 



A 

67 

15 

IF ( I CODE • EQ • 1 ) GO TO 1 



A 

68 


M5=UF/A5 



A 

69 


I 
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r 




7_0_ _ 

c 

OBTAIN 

MU5 ( US I NG RESULTS OF YOS ) AND RE5 FROM QDOT 

A 

71 

c 



A 

72 


CALL ODOT (0* ,0. «0. ) 

A 

73 

r 



A 

74 

r 

OBTAIN 

PREDICTED QT5 FROM QDOT IF INPUT QDO IS NOT 0 

A 

75 

c 



A 

76 


CALL QDOT ( 0« «0 • « 1 • ) 

A 

77 


PRINT 23 

A 

78 


PRINT 24 

A 

79 


PRINT 25, P5,RH05,T5,H5,S5R,Z5,GAME5,A5,U5,M5,RE5 

A 

80 


print 26 

A 

81 


PRINT 24 

A 

82 


PRINT 25, P5S,RH05S,T5S.H5S,S5SR,Z5S,GAME5S,A5S,U5S,M5S,RE5S 

A 

83 


print 27 

A 

84 


PRINT 28 

A 

85 


PRINT 29, PT5,RH0T5,T t 5,HT5,ST5R,ZT5, GAME T 5 , AT5 

A 

86 


IF (ODO 

• EQ. 0. ) GO TO 1 

A 

87 


PRINT 30 

A 

88 


PRINT 31 

A 

89 


PRINT 32, QT5CO,QT5HOS,QT5FR,QT5PT,QT5DE,QT5ZO,RN 

A 

90 


GO TO 1 


A 

91 

1 6 

STOP 


A 

92 

C 



A 

93 

1 7 

FORMAT 

(1 HI * A1 0// ) 

A 

94 

_1 8 

format 

(44H REAL-AIR DATA REDUCTION PROGRAM OF MILLER) 

A 

95 

1° 

FORMAT 

</53h all physical quantities in mks units- nasa sp-7012) 

A 

96 

?n 

FORMAT 

( /// 1 7H MEASURED INPUTS) 

A 

97 

2i 

format 

( /79H RUN PI VI PT QT R 

A 

9e 


I HOI 

RHOT TIMF) 

A 

99 

22 

FORMAT 

(BEIO.3 ) 

A 

100 

_23__ 

format 

(///24H FREE-STREAM CONDITIONS) 

_A_ 

JLOl 

24 

format 

C/l 07H P RHO T H S/R 

A 

102 


1 z 

GAME A V M NRF ) 

A 

1 03 

25 

format 

C 1 1F1 0 .3 ) 

A 

1 04 

26 

format 

(///39H STATIC CONDITIONS BEHIND NORMAL SHOCK) 

A 

1 05 

27 

FORMAT 

<///43H STAGNATION CONDITIONS BEHIND NORMAL SHOCK) 

A 

1 06 

28. 

FORMAT 

(/75H P RHO T H S/R 

A 

1 07 


1 z 

game a ) 

A_ 

1 08 

?Q 

FORMAT 

(8E1 0,3 ) 

A 

1 09 

30 

format 

C /// 44 H STAGNATION POINT HEAT TRANSFER PREDICTIONS) 

A 

1 10 

31 

format 

(/65H OT CO QTHOS QTFR QTPT QT DE Q 

A 

1 1 1 


1 TZ0 

RN ) 

A 

1 12 

32 

FORMAT 

(7F1 0.3 ) 

A 

1 1 3 
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END 

A 

1 14- 


SUBROUTINE PR0C2 ( PTsM , QT5M , XX ) 

B 

1 


COMMON /BLK 1 / RH05 . , S5R . T5 * A5 * Z5 , GAME5 , ZST AR5 * I CODE 

B 

2 


COMMON /BLK2/ RH0T5 . HT5 « ST5R . TT5 « AT5 * 2T5 . GAMET5 * ZS TART5 . PT5 

B 

3 


COMMON /BLK 3/ TWE • TAIJ# TW 4 QT5DE • SPH • TWALL • ETA • SAV 

B 

4 


COMMON /BLK4/ NVi IT iRiMMi I SP 1 NCO 1 MNi GAME5S 

B 

5 


COMMON /BLK5/ P5 . U5 . MB . MU5 . RE5 . Q T5FR » QT5H0S * QT5C0 . QT5Z0 , QT5PT 

B 

6 


COMMON /BLK6/ P5S iT5^i RHOSS t H5S iA5SiZ5S» UBS * MBS « MU5S * RE5S * S5SR 

B 

7 


COMMON /BLK7/ RN , TOLPT , TOLQT , TOLRHO , VOS . QDO . I TEST 

B 

8 


REAL M5* MU5f MSSfMUBS ,MUW#MUT5i NUtLAM 

B 

9 

c 


B 

10 

c 

FOR I TEST = 2 * XX=U5M AND USE DIRECT 

B 

I 1 

c 

FOR I TEST =3 « XX = P5M AND USE INVERSE 

B 

12 

c 


B 

13 


PT5=PT5M 

B 

14 


QT5=0T5M 

B 

15 

c 


B 

1 6 

c 

TWALL MUST BE 0 (TW=TwE ) OR BE 1 

B 

17 

r 


B 

18 

c 

FINITE VALUE OF ETA ( =CP*RHO*K ) MUST BE FURNISHED 

B 

19 


TW=TWE+1 • 77 2B-X-QTB* TWALL* SORT (TAU/ETA ) 

B 

20 


HW=1 .0046E+3*TW 

B 

21 


HT5=2.57745E+3*QT5*SQRT (RN/PT5 )*2.** (SPH/2. )+HW 

B 

22 


IF ( I TEST • EQ • 3 ) GO TO 1 

B 

23 

c 


. B 

24 

c 

FOR I TEST = 2 1 ESTIMATf INITIAL BO BY SETTING P5 = 0 

B 

25 

c 


B 

26 


P5 = 0. 

B 

27 


GO TO 2 

B 

28 

1 

P5M=YX 

B 

29 


P5=PBM 

B 

30 

2 

IF (SAV.EQ. 1 . ) GO TO 3 

B 

31 


CALL SEARCH (PT5 , RH0T5 , HT5 . ST5R . TT5 * AT5 . ZT5 . GAMET5 . ZSTART5 . ISP. 1 ) 

B 

32 


GO To 4 

B 

33 

3 

CALL SAVE ( PT5 . RH0T5 *HT5.ST5R.TT5.AT5.ZT5» GAMET5 . 3 ) 

B 

34 

r 


B 

35 

c 

IF TT5 .GE. 4500K, USE ZOBY RELATION TO OBTAIN HT5 

B 

36 

c 


B 

37 

4 

IF (TT5.GE.4500. ) GO TO 6 

B 

38 

C 


B 

39 

C 

IF TT5 .LT. 4500K * USE FAY-RIDDELL WITH MUT5 FROM YOS 

B 

40 

C 


B 

41 

5 

CALL QDOT ( QT5M * 1 • « 1 • > 

B 

42 

6 

IF ( I TEST • EQ • 3 ) GO TO 10 

B 

43 
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USM=XX 

R 

44 


U5 = URM 

B 

45 


H5=HT5-*5*U5**2 

B 

46 


IF CH5.GT.1 .E+5 ) GO TO 7 

B 

47 


PRINT 12 

R 

48 


ICODF=l 

B 

49 


GO TO 11 

B 

50 

7 

R5EST=PT5/( • 965*U5**? ) 

B 

51 


RH05=R5EST 

B 

52 


IF (FAV.EQ. 1*) GO TO 8 

B 

53 


CALL SEARCH < P5 , RHOS . H5 , S5R . T5 , A5 , ZB , GAMES , ZST ARB , I SP . 3 ) 

B 

54 


GO TO 9 

B 

55 

a 

CALL SAVE (P5,RH05,HS,S5R,T5,A5,Z5,GAME5,4) 

B 

56 

9 

CALL DIRECT (PT5M) 

B 

57 

c 


B 

58 

c 

BOCOM is COMPARISON op INITIAL BO TO UPGRADED B0(P5 ,NE. 0) 

B 

59 

c 


B 

60 


BOCOM=SQRT (PT5/ (PT5-P5 ) ) 

B 

61 

c 


R 

62 

c 

IF INITIAL BO NOT WITHIN 4 PERCENT OF UPGRADED BO, REPEAT 

B 

63 

c 


B 

64 


IF (BOCOM, LF". 1 .02 ) GO TO 11 

B 

65 


GO TO 5 

R 

66 

1 0 

CALL INVERSE (XX) 

B 

67 

1 1 

RETURN 

B 

68 

r 


B 

69 

12 

FORMAT (46H HI IS LESS THAN MINIMUM VALUE ON TAPE , I TEST = 2 ) 

B 

70 


END 

B 

71- 


SUBROUTINE PR0C4 ( P5M , U5M , QT5M ) 

C 

1 


COMMON /BLK1/ RH05,H=;,S5R,T5*A5,Z5,GAME5,ZSTAR5»I CODE 

C 

2 


COMMON /BLK2/ RH0T5 , HT5 , ST5R , TT5 , AT5 , ZT5 , GAMET5 , ZS TART5 # PT5 

C 

3 


COMMON /BLK3/ TWE , T Al Jl TW , QT5DE , SPH , TWALL , ETA , SAV 

C 

4 


COMMON /BLK4/ NV t I T # R • MM • I SP * NCO * MN « GAME5S 

c 

5 


COMMON /BLK5/ P5 , U5 , M5 , MU5 , RE5 , QT5FR , QT5H0S , QT5C0 . QT5Z0 , QT5PT 

c 

6 


COMMON /BLK6/ P5S , T5S , RH05S , H5S , A5S , Z5S , U5S , MBS , MU5S , RE5S , S5SR 

c 

7 


COMMON /BLI< 7/ RN , TOLPT , TOLQT , TOLRHO , YOS , QDO , I TEST 

c 

8 


REAL M5,MU5,MBS,MU5S,MUW,MUT5, NU , LA M 

c 

9 


NN=n 

c 

1 0 


OTE=OTEM 

c 

1 1 


P5=PqM 

c 

12 


U5=U^M 

c 

1 3 

C 


c 

14 

r 

TWALL must BE 0 (TW=TI.(E ) OR BE 1 

c 

15 

C 


c 

1 6 
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c 

FINITE VALUE OF ETA < =CP*RHO*K ) MUST BE FURNISHED 

TW = TWE+1 .7725*QT5*TWALL*SQRT (TAU/FTA ) 

c 

C 

1_7__ 

18 

■ 

H W = 1 • 0046E+3*TW 

R5EST=2 *684 1 E + 7*RN-*0t5M**2/U5M**6 

RHOUD=2 .05*R5EST 

C 

c 

c 

21 


RH0LOW=.20*R5FST 

c 

22 

r 


c 

23 

C 

USFR CAUTIONED ON THpSE LIMITS OF RH05(RH01) 

c 

24 

C 

THESF limits SHOULD pE valid for nearly all cases 

c 

25 

c 

TO MINIMIZE COMPUTER TIME, SHOULD REFINE LIMITS 

c 

26 

c 

rOR M5 GREATER THAN 10, LFT RHOLOW BE .85*R5EST 

c 

27 

r 

AND RHOUP BF 1.05*R5EST 

c 

28 

C 


C 

29 


RH05=1 • 45*R 5EST 

C 

30 

I 

DFLRhO= (RHOuP-RHOLOUI )/ 4 . 

c 

31 


IF (NN.EQ.O) "go TO 2 

c 

32 


RH05=RH0UP-DELRH0 

c 

33 

2 

IF (bAV.EQ. 1 . ) GO TO 3 

c 

34 


CALL SEARCH ( P5 « RH05 , H5 , S5R » T5 . A5 , Z5 . GAME5 , ZS T AR5 * I SP , 2 ) 

c 

35 


GO TO 4 

c 

36 

3 

CALL SAVE (P5,RH05iH=;,S5R,T5.A5.Z5.GAME5.2 ) 

c 

37 

4 

PT5=.965*RH05*U5**2 

c 

38 


PT5M =PT5 

c 

39 


CALL DIRECT (PT5M) 

c 

40 


IF fTTB.LT.4B00.) GO TO 6 

C 

41 


QTBZO = 3 . 8798E — 4*SQRT ( PTB/RN ) * ( HT5-HW ) /2 . ** (SPH/2. ) 

c 

42 

r 


c 

43 

c 

TOLOT REPRESENTS TOLERANCE OF ITERATION ON QT5 

c 

44 

r. 


C 

45 


IF ( ABS ( I .-QT5M/QT5Z0 ) . LE .TOLQT ) GO TO 8 

c 

46 


IF (MN.E0.15) GO TO 7 

C 

47 


NN = NN+ 1 

c 

48 


IF (0T5Z0.LT. QT5M ) GO TO 5 

c 

49 


RH05=RH05-DELRH0 

c_ 

50 


GO T n 2 

c 

__51 

5 

RHOLOW =RH05 

c 

52 


RH0UD=RH05+DELRH0 

c 

53 


GO To l 

c 

54 

6 

CALL QDOT ( QT5M , 1 • , 1 * ) 

c 

55 

c 


c 

56 

c 

TOLQT REPRESENTS TOLFRANCE OF ITERATION ON QT5 

c 

57 

c 


C 

58 


IF ( ABS ( 1 .-QT5M/QT5FR) .LE. TOLQT ) GO TO 8 

c 

59 


IF (NN.EQ.15) GO TO 7 

c 

60 
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NN=NM+ 1 

C 

(SI 


IF (OT5FR.LT.QT5M) GO TO 5 

c 

62 


RH05=RH05— DELRHO 

c 

63 


GO TO 2 

c 

64 

7 

PRINT 9 

c 

65 


PRINT 10 

c 

66 


IC0DP=1 

c 

67 

8 

RETURN 

c 

68 

0 


c 

69 

9 

FORMAT (50H ITERATIONS ON RH05 EXCEED LIMIT OF 15 FOR ITEST=4) 

c 

70 

1 0 

FORMAT ( 2 1 H CHECK LIMITS ON RH05 ) 

c 

71 


FND 

c 

72- 


SUBROUTINE PR0C5 ( RH05M ♦ P5M * Y Y ) 

D 

1 


COMMON /BLK 1 / RH05 » H5 i S5R * T5 i A5 * Z5 . GAME5 * ZST AR5 i I CODE 

D 

2 


COMMON /BLK2/ RH0T5 * HT5i ST5R » TT5 » AT5 » ZT5 »GAMET5iZSTART5*PT5 

D 

3 


COMMON /BLK 3/ TWE*TAIliTW* QT5DE iSPHt TWALL 1 ETA 1 SAV 

D 

4 


COMMON /BLK 4/ NV ♦ I T % R * MM * I SP ♦ NCO * MN , GA ME5S 

D 

5 


COMMON /BLK 5/ P5,U5iM5, MU5 • RE5 < QT5FR , QT5H0S * QT5C0 t QT5Z0 « QT5PT 

D 

6 


COMMON /BLK6/ P5S * T5S ♦ RH05S * H5S« A5S ♦ Z5S n U5S *M5SiMU5S» RE5S « S5SR 

D 

7 


COMMON /BLK7/ RN , TOLPT , TOLQT , TOLRHO . YOS , QDO » I TEST • 

D 

8 


REAL M5 « MU5 * M5S • MU5S ,MUW * MUT 5 * NU * L AM 

D 

9 

C 


D 

1 0 

c 

FOR I TEST = 5 • YY = U5M AND USE DIRECT 

D 

1 1 

c 

FOR I TEST = 6 . Y Y = PT5M AND USE DIRECT 

D 

12 

c 


D 

1 3 


RH05=RH05M 

D 

14 


P5=PFM 

D 

15 


IF (ITEST.EQ.5) GO TO 1 

D 

1 6 


PTBMrrYY 

D 

1 7 


ptb=ptem 

_D_ 

18 


GO TO 2 

D 

19 

1 

U 5 M = Y Y 

D 

20 


UB =U C 'M 

D 

21 

2 

I F ( ^AV • EO • 1 .) GO TO 3 

D 

22 


CALL SEARCH ( P5 .RH05 , H5 , S5R ♦ T5 i A5 , Z5 . GAME5 , ZST AR5 . I SP , 2 ) 

D 

23 


GO TO 4 

D 

24 

3 

CALL SAVE <P5 .RH05.H5. S5R.T5 . A5» Z5.GAME5.2 ) 

D 

25 

4 

IF ( I TFST • EG • 5 ) GO TO 5 

D_ 

26 

r 


0_ 

27 

c 

FOR I TEST = 6 » MUST ESTIMATE U5 EOR DIRECT 

D 

28 

c 

UPGRADING OF U5 PERFORMED IN DIRECT 

D 

29 

r- 


D_ 

30 


U5=S0RT ( PT5/ ( ,965*RH05 ) ) 

D 

31 


GO TO 6 

D 

32 
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c 


D 

33 

c 

FOR I TEST = 5 , MUST ESTIMATE PT5M FOR USE IN DIRECT 

D 

34 

c. 


D 

35 

E 

PT5= # 965*RH05*U5**2 

D 

36 


PT5M=PT5 

D 

37 

6 

CALL DIRFCT ( PT5M ) 

D 

38 


RFTURN 

D 

39 


FND 

D 

40 — 


SUBROUTINE PR0C7 ( PTSM , RH0T5M , YY ) 

E 

1 


COMMON /BLK 1 / RH05,H^,S5R,T5.A5,Z5,GAME5,ZSTAR5, ICODE 

E 

2 


COMMON /BLK2/ RH0T5 , HT5 , ST5R , TT5 , AT5 , ZT5 , GAMET5 , ZSTART5 , PT5 

. E 

3 


COMMON /BLK3/ TWE , T AU « TW , QT5DE , SPH , TWALL , ETA , S A V 

E 

4 


COMMON /BLK4/ NV • I T * P < MM , I SP * NCO , MN , GAME5S 

E 

5 


COMMON /BL<5/ P5,U5,M5,MU5,RE5,QT5FR,QT5H0S,QT5C0,QT5Z0.QT5PT 

E 

6 


COMMON /BLK6/ P5S , T5S , RH05S , H5S , A5S , Z5S , U5S , M5S , MU5S , RE5S , S5SR 

E 

7 


COMMON /BLK7/ RN , TOLPT , TOLQT , TOLRHO , VOS , QDO , I TEST 

E 

8 


REAL M5,MU5,M5S,MU5S,MUW,MUT5,NU,LAM 

E 

9 

C 


E 

10 

C 

FOR I TEST = 7 % YY = RH05m AND USE INVERSE 

E 

1 1 

C 

FOR I TEST = B « YY = P5M AND USE INVERSE 

E 

32 

C 

FOR I TEST =9 i YY = U5M AND USE DIRECT 

E 

13 

C 


E 

14 


PT5=PT5M 

E 

15 


RH0TE=RH0T5M 

E 

1 6 


IF (SAV.EQ. 1.) GO TO 1 

E 

17 


CALL SEARCH ( PT5 , RH0T5 , HT5 * S T5R . TT5 * AT5 , ZT5 * GAMET5 , ZSTART5 i ISP*2) 

E 

18 


GO TO 2 

E 

19 

1 

CALL SAVE (PT5.RH0T5«HT5«ST5R.TT5.AT5.ZT5.GAMET5 .2) 

E 

20 

c 


E 

21 

c 

at this point, know stagnation conditions exactly 

E 

22 

c 


E 

23 

2 

IF ( r TEST .EQ. 9 ) GO To 3 

E 

24 


CALL INVERSE (YY) 

E 

25 


GO TO 7 

E 

26 

3 

U5M=YY 

E 

27 


U5=URM 

E 

28 


H5=Ht5-.5*U5**2 

E 

29 


IF (H5.GT, 1 .E + 5 ) GO TO 4 

E 

30 


PRINT 8 

E 

31 


ICOD^=l 

E 

32 


GO TO 7 

E 

33 

4 

R5EST=PT5/( .965*U5**?> 

E 

34 


RH05=R5EST 

E 

35 


IF (SAV.EQ. 1 . ) GO TO 5 

E 

36 
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05 



CALL SEARCH (PS • RH05 • H5 • S5R • T5 • AS • Z5 • GAMES . ZSTAR5 • I SP • 3 ) 

F 

37 


GO TO 6 

E 

38 __ 

5 

CALL SAVE (P5,RH05,HB,S5R,T5,A5,Z5,GAME5,4) 

E 

39 

6 

CALL DIRECT ( PT5M ) 

E 

40 

7 

return 

E 

41 

C 


E 

42 

8 

format (50H hi IS LESS THAN MINIMUM VALUE ON TAPE FOR ITEST=9> 

E 

43 


END 

E 

44- 


SUBROUTINE DIRECT (PT=5M) 

F 

1 


DIMENSION X<4)* Y(4*Q*150)* Z(9)* U(4)* V(4)» W(4)* NP ( 4 ) 

F 

2 


DIMENSION T ABP ( 3 ) « T ABH ( 3 ) 

F 

3 


COMMON /BLK 1 / RH05.Hf5»S5R,T5*A5»Z5*GAME5*ZSTAR5* I CODE 

F 

4 


COMMON /B.-K2/ RH0T5 , HT5» ST5R »TT5»AT5»ZT5,GAMET5*ZSTART5,PT5 

F 

5 


COMMON /BLK3/ TWE»TAu*TW*QT5DE»SPH«TWALL*ETA«SAV 

F 

6 


COMMON /BLK4/ NVi IT iRtMUi ISPi NCO iMNiGAME5S 

F 

7 


COMMON /BLK5/ P5 » U5 , M5 * MU5 , RE5 . QT5FR , QT5H0S , QT5C0 * QT5Z0 , QT5PT 

F 

8 


COMMON /BLK6/ P5S»T5SiRH05S,H5S«A5S*Z5S*U5S*M5S.MU5S.RE5S*S5SR 

F 

9 


COMMON /BLK7/ RN, TOLPT.TOLQT , TOLRHO i YOS.QDO * I TEST 

F 

1 0 


COMMON I COUNT * I MET ( 2 ) » NP» ABARi ME »MF 

F 

1 1 


REAL M5i MU 5 iM5Si MU5S « MUW , MUT5 « NU » LAM 

F 

12 


IMET ( 1 ) = IMET <2 ) =0 

F 

1 3 

C 


F 

1 4 

c 

DIRECT PERFORMS NORMAl SHOCK CROSSING* PRE-TO-POST 

F 

15 

C 


F 

1 6 

1 

BSNS=RH05*U5 

F 

1 7 


CSNS=P5+BSNS*U5 

F 

18 


DSNS=H5+.B*U5**2 

F 

1 9 


HT5=H5+*5*U5**2 

F 

20 


CALL SAVE (PT5,RH0T5,HT5,ST5R,TT5.AT5,ZT5,GAMET5,3) 

F 

21 


RH05S=* 955*RH0T5 

F 

22 

2 

U5S=RSNS/RH05S 

F 

23 


PSS=CSNS-BSNS*U5S 

F 

24 


H5S=nSNS-.5*U5S**2 

F 

25 


IF (SAV.FO. 1 . ) GO TO 3 

F 

26 


CALL SEARCH ( P5S , RNE'iU H5S , S5SR , T5S « ASS , Z5S , GAME5S . ZST AR5S , ISP, 1 ) 

F 

27 


GO TO 4 

F 

__28_ 

3 

CALL SAVE (P5S, RNEW, H5S, S5SR ,T5S * A5S , Z5S , GAME5S, 3 ) 

F 

29 

4 

IE ( ABS ( 1 • -RH05S/RNE"* ) • LE • • 00 1 ) GO TO 5 

F 

30 


RH05 S-RNEW 

F_ 

31 


GO TO 2 

F 

32 

5 

RH05F = RNEW 

F 

33 


M5S=uSS/A5S 

F 

34 


HT5R=HT5/R 

F 

35 


IF (SAV.EQ. 1 .) GO TO 6 

F 

36 
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XX = S e ?SR 

F 

37 


Z (4 ) =AL0G1 0 (HT5R) 

F 

38 


CALL SLOW (XX.Z.4.3, TT.NV.NERR.Y.X ) 

F 

39 


PT5=(10#**Z(3) >*1 .0I325E+5 

F 

40 


GO TO 8 

F 

41 

6 

PT5 = P5S*('l . + ( (GAME5S-1 . ) /?_ . ) *M5S**? ) ** ( GAME5S/ ( GAME5S- 1 . ) ) 

F 

42 


TABP ( 1 ) = ,95*PT5 

F 

43 


TABP (2 ) =PT 5 

F 

44 


TABP(3 ) =1 • 05*PT5 

F 

45 


DO 7 1=1,3 

F 

46 


PT5=TABP ( I ) 

F 

47 


CALL SAVE (PT5,RHOT5,HT5A,S5SR,TT55.AT5,ZT5,GAMET5» 1 ) 

F 

48 


T ABH ( I ) = HT5 A 

F 

49 

7 

CONT JNUF 

F 

50 


CALL FTLUP ( HT5 , PT5 , ? . 3 * TABH * T ABP ) 

F 

51 


CALL SAVE <PT5,RH0T5,HT5,S5SR,TT5,AT5.ZT5,GAMET5.1 ) 

F 

52 

8 

IF ( I TEST. EO. 4. OR. I TPST.EQ.5 ) GO TO 13 

F 

53 

C 


F 

54 

c 

TOLPT REPRESENTS TOLFRANCE OF ITERATION ON PT5 

F 

55 

c 


F 

56 


IF ( ABS ( I .-PT5M/PT5 ) ,LE. TOLPT ) GO TO 13 

F 

57 


IF ( ITEST.EQ.2.0R. ITFST.EQ.9 ) GO TO 10 

F 

58 


IF ( 1 TEST • EQ • 1 ) GO TO 9 

F 

59 


U 5=tic ? *SQRT (PT5M/PT5 ) 

F 

60 


GO TO 1 

F 

61 

Q 

H5=H^*PT5/PT5M 

F 

62 


GO TO 11 

F__ 

63 

in 

P5=P5*PT5M/PT5 

F 

64 

1 1 

IF ( c; A V • EQ • U) GO TO 12 

F 

65 


CALL SEARCH (P5.RH05.H5, S5R ,T5,A5,Z5, GAME 5 , ZST AR5 * I SP , 1 1 

F 

66 


GO TO 1 

F 

67 

12 

CALL SAVE (P5,RH05,H=;.S5R,T5«A5»Z5,GAME5.3) 

F 

68 


GO TO 1 

F 

69 

1 3 

IF (FAV.EO. 1 . ) GO TO 14 

F 

70 


CALL SEARCH (PT5 , RH0T5 , HT5 , ST5R * TT5 , AT5 . ZT5 , GAMET5 , ZSTART5 , ISP. 1 ) 

F 

71 


GO TO 15 

_F_ 

72 

14 

CALL SAVE (PT5.RH0TS,HT5.ST5R.TT5.AT5.ZT5,GAMET5.3) 

F 

73 

1 =5 

RFTURN 

F 

74 


END 

F 

75- 


SUBROUTINE INVERSE <7Z) 

G 

1 


DIMENSION X ( 4 ) . Y<4. 9.150), Z(9), U(4). V(4), W<4)» NP(4) 

G 

2 


COMMON /BLK 1 / RH05 , H 5 , S5R , T 5 , A5 , Z5 , GAME5 , ZST AR5 , I CODE 

G 

3 


COMMON /BLK2/ RH0T5,hT5,ST5R.TT5,AT5,ZT5,GAMET5,ZSTART5,PT5 

G 

4 


COMMON /BLK3/ TWE t TAlJ* TW * QT5DE ♦ SPH , T WALL t ETA , SAV 

G 

5 


-a 
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COMMON /BLK4/ NV » IT 1 R 1 HM 1 ISP* NCO ,MN,GAME5S 

G 

6 


COMMON /BLK5/ P5 • U5 , M5 , MU5 , RE5 , QT5FR , QT5H0S * QT5C0 , QT5Z0 , QT5PT 

G 

7 


COMMON /BLK6/ P5S,T5E,RH05S,H5S, A5S,Z5S,U5S,M5S,MU5S,RE5S,S5SR 

G 

8 


COMMON /BLK7/ RN,TOLPT , TOlQT , TOLRHO , YOS , QDO , ITEST 

G 

9 


COMMON I COUNT , I MET ( 2 ) , NP , ABAR , ME , MF 

G 

1 0 


REAL M5«MU5*MES*MU5S ,MUW« MUTSiNU • LAM 

G 

1 1 


IMFT ( 1 ) = I MET (2 ) =0 

G 

12 

c 


G 

13 

c 

INVERSE PERFORMS NORMAL SHOCK CROSSING, POST-TO-PRE 

G 

14 

Q 


G 

15 


XX=ST5R 

G 

1 6 __ 


MM = 0 

G 

17 ... 


PLOW=.850*PT5 

G 

16 


PUP= »970*PT5 

G 

19 

c 


G 

20 

r 

USER CAUTIONED ON THfSE LIMITS OF P5SIP2) 

G 

21 

c 

IF MF LESS THAN 3, MAY HAVE TO LOWER PLOW 

G 

22 

c 

TO MINIMIZE COMPUTER TIME, SHOULD ALSO LOWER PUP 

G 

23 

c 


G 

24 


P5S=PUP 

G 

25 

i 

DELP= ( PUP-PLOW ) /4 . 

G 

26 


IF (MM.EO.O ) GO TO 2 

G 

27 


P5S=PUP-DELP 

G 

_28 

__2_ 

IF (FAV.EQ. 1 . ) GO TO 3 

G 

29 


Z (3 )=AL0G1 0 (P5S/1 .01325E+5) 

G 

30 


CALL SLOW ( XX ,Z , 3,4 , IT ,NV ,NERR, Y ,X > 

G 

31 


CALL SLOW (XX, Z, 3,2, IT ,NV, NERR , Y , X ) 

G 

32 


H5S = ( 1 0.**Z ( 4 ) >*R 

G 

33 


RH05F= ( 1 J » **Z ( 2 > )*1 ,P9 14889 

G 

34 


GO TO 4 

G 

35 

3~ 

CALL SAVE (P5S,RH05S,H5S,ST5R,T5S,A5S«Z5S,GAME5S, 1 ) 

G 

36 

4 

U5S= FORT ( 2 , * ( HT5-H5S ) ) 

G 

37 


IF (MM. EG. 38) GO TO 13 

G 

38 


MM=MM+1 

G 

39 


IF ( ITEST. EQ. 3. OR. ITFST.EQ.8 ) GO TO 9 

G 

40 


RH05M=ZZ 

G 

41 


RH05=RH05M 

G 

42 


U5=RH05S*U5S/RH05 

G 

43 


P5=PeS+RH05S*U5S**2-RH05*U5**2 

G 

44 


IF ( P5 . LT • • 1 ) GO TO fl 

G 

45 


he = hte- . 5 * 115**2 

G 

46 


IF (H5.LT. 1 .E + 5 ) GO TO 8 

G 

47 


IF (CAV.FQ. 1 . ) GO TO 6 

G 

48 


CALL SEARCH (P5 « RHOE ,H5 , S5R ,T5,A5,Z5, GAME 5 , ZST AR5 , I SP , 1 ) 

G 

49 
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GO TO 7 

G 

50 

6 

CALL SAVE (P5,RH0E,H^.S5R,T5.A5.Z5,GAME5»3) 

G 

51 

C 


G 

52 

C 

TOLRHO REPRESENTS TOLERANCE OF ITERATION ON RH05 

G 

53 

c 


G 

54 

7 

IF ( ABS ( 1 .-RH05/RH0E ) .LE. TOLRHO ) GO TO 10 

G 

55 


IF fPH0E.LT.RH05) GO TO 8 

G 

56 


P5S=P5S“DFLP 

G 

57 


GO TO 2 

G 

58 

B 

PL0W=P5S 

G 

59 


PUP=P5S+DFLP 

G 

60 


GO TO 1 

G 

61 

9 

P5M=7Z 

G 

62 


P5=PFM 

G 

63 


U5= (P5S-P5 )/ <RH05S*UFS )+U5S 

G 

64 


RH05=RH05S*U5S/U5 

G 

65 


H5=Ht 5— . i*U5**2 

G 

66 


IF (H5.GT. 1 .E+5 ) GO TO 5 

G 

67 


PES=P5S— DELP 

G 

68 


GO To 2 

G 

69 

10 

IF (FAV.EQ. 1 . ) GO TO 11 

G 

70 


CALL SEARCH ( P5S » RHO5S » H5S * S5SR »T5S»A5S*Z5S*GAME5S» ZSTAR5S * ISP* 1 ) 

G 

71 


GO TO 12 

G_ 

72 

1 1 

CALL SAVE (P5S,RH05S,H5S*S5SR.T5S,A5S.Z5S,GAME5S.3 ) 

.G 

73 

12 

M5S=i_)5S/A5S 

G 

74 


GO TO 14 

G 

75 

13 

PRINT 15 

G 

76 


IC0DP=1 

G 

77 

14 

return 

G 

78 

C 


G 

79 

1=5 

FORMAT (52H ITERATIONS ON P2 EXCEED LIMIT - REFINE LIMITS ON P2 ) 

G 

80 


END 

G 

81- 


SUBROUTINE QDOT ( QTBM , B I T , D I T ) 

H 

1 


DIMENSION TABP(7), T ABT ( 25 ) * TABNU(175) 

H 

2 


DIMENSION X ( 4 ) * Y <4* <3. 1501* Z(9), U<4), V(4)« W<4)* NPC4) 

H 

3. 


DIMENSION TAPY(4), T ABTY (13 1* T ABNU Y ( 52 ) 

H 

4 


COMMON /BLK 1 / RH05,H5*S5R,T5,A5.Z5,GAME5.ZSTAR5. ICODE 

H 

5 


COMMON /BLK2/ RH0T5 * HT5 * ST5R • TT5 i AT5 ♦ ZT5 * GAMET5 * ZSTART5* PT5 

H 

6 


COMMON /BLK3/ TWE . TAU * TW ♦ QT5DE . SPH , TWALL * ETA , S AV 

H 

7 


COMMON /BLK4/ NV. I T , R * MM , I SP * NCO , MN * GAME5S 

H 

8 


COMMON /BLK5/ P5 , U5 . ME , MU5 . RE5 * QT5FR * 0T5H0S • QT5C0 * QT5Z0 * QT5PT 

H 

9 


COMMON /BLK6/ P5S . T5F , RH05S , H5S « A5S « Z5S , U5S * M5S « MU5S . RE5S * S5SR 

H 

10 


COMMON /BLK7/ RN , TOLPT , TOLQT , TOLRHO * YOS , QDO « I TEST 

H 

1 1 


REAL M5,MU5,M5S*MUBS,MUW*MUT5*NU.LAM 

H 

12 


APPENDIX D - Continued 



C7I 

O 



r 


H 

13 

C 

TABL^ OF VISCOSITY FROM HANSEN ( NASA TR R-50 ) 

H 

14 

c 


H 

15 


DATA TABP/1 .01 325E+7, 1 .01325E+6. 1 .01325E+5, 1 .01325E+4, 1 .01325E+3, 1 

H 

16 


1 . 0 1 325E+ 3 « 1 .01 325E+ 1 / 

H 

17 


DATA TABT/3000. , 3500 . , 4000 . * 4500 . , 500C . , 5500 . , 6000 . .6500. . 7000. ,75 

H 

18 


100. ,8000.48500.4 9000.* 9500. % 1 0000. , 1 0500 • i 1 1 000. ♦ 1 1 500 • , 1 2000 • ♦ 1 25 

H 

19 


200 . , 13000. . 13500. , 14H00. . 14500. . 1 5000./ 

H 

20 


DATA TABNU/2*1 . 0. 1. 0 03 . l.olO. 1. 02E, 1 .036, 1 .050 ,1.072,1.089,1.112,1 

H 

21 


1.143,1.185,1 .238, 1 ,2q8 , 1 . 36 1 , 1 . 4 1 8 , 1 . 467 , 1 .509,1 .549,1 .577, 1.581,1 

H 

22 


2.594 , 1 .599, 1.601 ,1 .604 ,1.0,1.001,1.008,1. 022 ,1.036,1.052,1 .067, 1 .0 

H 

23 


390,1.124,1.175,1.238,1.307,1 .368 ,1.418,1.458,1.496,1.501,1.511,1.5 

H 

24 


420,1.516,1.508,1.4 92,1 .468,1 .415,1 .387 ,1., 1.003, 1.016, 1.0 29, 1,043* 

H 

25 


51 . 06 n , 1.090,1.139,1.?08,1.283,1 .342, 1 .386, 1 . 425 , 1 . 438 , 1 . 445 , 1 .448 , 

H 

26 


61 .442, 1 *424 , 1 .394 ,1.342,1 .274 , 1 . 1 87 , 1 . 082 , . 94 , . 828 , 1 ♦ , 1 . 006 , 1 . 02 * 1 

H 

27 


7.033, 1.051,1 .086, 1.148,1 .229,1 .294, 1 .332 ,1.371,1 .386, 1 .396, 1 .393, 1 

H 

28 


8.375, 1 • 335 , 1 .267 ,1 .168,1.0404.881 ,.7114.5474.408, .268, .21241. .1.01 

H 

29 


9,1.022.1.038,1 .074 , 1 . 1 46. 1 .228, 1 .276. 1.317,1 .337, 1 .347. 1 .343, 1.314 

H 

30 


S , 1 .2^1 , 1 . 1 43 4 .983 4 . 732 4.571 4 .387 4 • 249 4 • 1 58 , • 1 00 4 • 067 4 .042 . *0 1 6 * 1 • * 

H 

31 


Sl .01 , 1 . 024. 1.055,1.128,1.209,1.257,1 .286,1 .303, 1 . 307. 1 . 28 , 1 .207 , 1 . 

H 

32 


S068 4 .853 4 .595 4 .361 4 • ?CO 4 • 1 08 4 .063 4.036 4 .024 4 .0184 .0154.0134.01241 • 

H 

33 


S, 1.011,1. 032, 1.096, 1.181, 1.2 27, 1.256,1. 271,1 .264 ,1.210, 1.072, .826. 

H 

34 


S. 517, .261 4. 118, .0 55 4 ,029, .018,. 012,. 0 09,. 008, .007, .007 4.0 03, .008/ 

H 

35 

r 


H 

36 

C 

TABLF OF VISCOSITY FROM YOS ( AVCO RAD-TM-63-7 > 

H 

37 

r 


H 

38 


data TAPY/1 .01 325E+5.3. 03975E+5 , 1.01 325E+6 , 3 . 0 3975E+6/ 

H 

39 


DATA TABTY/1000. ,2000 . , 3000. .4000. ,5000. ,6000. ,7000. ,8000 . .9000. . 1 

H 

40 


1 0000 # f 1 2000 • 4 1 4000. ,16000./ 

H 

41 


data TABNUY/.4 1 8E— 4 , .648E-4, .858E-4, 1 .08E-4, 1 .30E-4.1 .54E-4, 1 .86E- 

H 

42 


1 4 ,2.21 E-4 , 2.46E-4 ,2 . 63E-4 , 2 . 53E-4 , 1 . 77E-4 , .96E-4 , . 418E-4 . .648E-4 , . 

H 

43 


2857E-4 , 1 • 07E-4 , 1 .30E-4 , 1 .52E-4, 1 .80E-4, 2. 1 4 E-4 , 2 .45E-4 , 2 . 66E-4 , 2 .8 

H 

44 


35E-4 , 2 . 34E-4 , 1 .53E-4 , .41 8E-4 , • 64 8E-4 , • 857E-4 , 1 . 07E-4 , 1 .30 E-4 , 1 • 51 E 

H 

45 


4-4,1 .76E-4. 2.06E-4,2.4E-4,2.67E-4 .3.00E-4 .2.82E-4 , 2.24E-4 , .418E-4, 

H 

46 


5.648=--4, • 856E-4 , 1 .06F-4, 1 .27E-4, 1 .50E-4, 1 . 73E-4 , 2 . OOE-4 , 2 . 32E-4 ♦ 2 . 

H 

47 


663 E-4 4 3 • 06E — 4 , 3 • 1 OE- 4 , 2 • 66E- 4/ 

H 

48 


MM = 0 

H 

49 

c 


__H_ 

50 

c 

QDOT CONSISTS OF 3 SECTIONS 

H 

51 

" c 


H 

52 


l<=- (niT.EQ. 1 .) GO TO 5 

H 

53 

r 


H 

54 

C 

(1) PREDICTING MU5 < I NG YOS RESULTS) AND RE5 

H 

55 

C 


H 

56 


APPENDIX D - Continued 



IF (tF.LF. 1 500. ) GO TO 1 

H 

57 


CALL DISCOT ( T5 . P5 . TAP TY . T ABNUY . T AP Y , 1 1 , 52 , 4 , MU5 ) 

H 

58 


GO T n P 

H 

59 

1 

MU5=1 . 462 E- 6* SORT ( T5 ) / ( 1 . + 1 12./T5) 

H 

60 

2 

RE5 = PH05*U5/MIJ5 

H 

61 


IF (T5S.LF. 1500. ) GO TO 3 

H 

62 


CALL CISCOT (T5S.P5S.TABTY.TABNUY.TAPY. 1 1 .52.4. MU5S) 

H 

63 


GO TO 4 

H 

64 

3 

MU5S=1 . 4 62 E- 6* SORT (TBS )/( 1 . + 1 12 ./T5S ) 

H 

65 

4 

RF5S=RH05*U5/MU5S 

H 

66 


GO TO 20 

H 

67 

n 


H 

68 

c 

(2) PREDICTING 0T5 WITH HANSEN OR YOS RESULTS FOR MUT5 

H 

69 

r 


H 

70 

5 

PR=.7l 

H 

71 


NN = 0 

H 

72 


TW=TWF 

H 

73 

6 

HW=1 . 0046F+3*TW 

H 

74 


RHOW=3 • 48398E— 3*PT5/TW 

H 

75 


MUW= 1 . 462E-6*SQRT (TW )/ < 1 . + 1 1 2./TW ) 

H 

76 


BO= (1 ./RN )*SQRT (2 .* (PT5-P5 1/RHOT5 ) 

H 

77 


IF (RIT.FQ. 1 . ) GO TO 16 

H 

78 

C 


H 

79 

C 

IF 000=0. VARIOUS QTR WILL NOT BE COMPUTED 

H 

80 

r 


H 

81 


IF (OOO.FQ.O.) GO TO 20 

H 

82 


IF (TT5.GT. 1500. ) GO TO 7 

H 

83 


MUT5= 1 .462E— 6*SQRT (Tt5 1/(1. +112. /TT5) 

H 

84 


C-0 TO P 

H 

85 

r 


H 

86 

c 

IF YOS=0 . USE HANSENr RESULTS FOR MUT5 

H 

87 

c 

IF YOS= 1 . USE YOS RESULTS FOR MUT5 

H 

88 

r 


H 

89 

T 

IF (YOS.EQ.O. ) GO TO 8 

H 

90 


CALL DISCOT ( TT5 , PT5 , T ABTY . T ABNUY . TAPY , 1 1 . 52 , 4 , MUT5 ) 

H 

91 


GO TO 9 

H 

92 

8 

CALL DISCOT (TT5.PT5 ,TABT, TABNU. TABP. 1 1 . 1 75.7.NU 1 

H 

93 


MUT5= 1 .462E-6*SQRT ( TT5 ) / ( 1 .+ 1 1 2 ./TT5 ) 

H 

94 


MUT5=MUT5*NU 

H 

95 

9 

QT5H0S= 1 . 1 672* (SORT ( BO*RHOW*MUW*TW** . 4 ) )* (U5**l .69 )* ( 1 .-HW/HT51/2 . 

H 

96 


1 ** (SPH/2. ) 

H 

97 


QT5FR= .760* ( ( RHOW*MUW ) ** . 1 )* ( (RH0T5*MUT5 )**.4 ) * ( HT5-HW ) *SQRT (BO )/ ( 

H 

98 


1 <PR**0.6)*2.**(SPH/2. ) ) 

H 

99 


IF (U5.LE.8. 84E+3 ) 10.11 

H 

100 
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i 


1 0 

FPS=n . 0 

H 

1 01 


GO TO 12 

H 

102 

1 1 

FPS=1 . 

H 

103 

1 2 

QT5CO=.7670* ( (RHOW*Ml)W )#*.07)*< <RH0T5*MUT5 ) ** . 43 ) * ( HT5-HW )*SQRT (BO 

H 

104 


1 )* ( I . +. 075*FPS* ( HT5/1 • 9686E+ 7-2 . >**2 )/ ( ( PR**0 • 6 1 *2 . ** ( SPH/2 . ) ) 

H 

105 


QT5ZO=3.8798E-4*SQRT ( PT5/RN ) * ( HT5-HW ) /2 . ** (SPH/2 . ) 

H 

106 


IF (lj5.LF*9.906E + 3> 13*14 

H 

107 

1 3 

LAM= 1 • 

H 

108 


GO TO 15 

H 

109 

1 4 

l_AM=u5/9.906F+3 

H 

1 10 

1 5 

QT5PT =0.9000* ( (RH0W*MUW )**.07 )* ( <RH0T5*MUT5 )**.43 ) * (HT5-HW ) *SQRT ( B 

H 

1 1 1 


1 0 )*LAM/ ( (PR**0 .75 )*2.** ( SPH/2. ) ) 

H 

1 1 2 


QTFiDP = • 62* ( (RHOW*MUW ) **• 1 624 )* ( ( RHOT5*MUT5 ) ** • 3376 ) * < HT5-HW ) * SORT ( 

H 

113 


1 BO )/(PR*2.**( SPH/2. ) ) 

H 

114 


IF (TWALL •FQ.O. ) GO TO 20 

H 

1 15 


IF ( \iN • E O # 2 ) GO TO 2 n 

H 

1 16 


NN=NM+ 1 

H 

117 


QTB=OT5ZO 

H 

1 18 


TW=TwE+l «7725*QT5*TWALL*SQRT (TAU/ETA ) 

H 

1 19 


GO To 6 

H 

120 

c 


H 

121 

C 

(3) PREDICTING HT5 FROM FAY-R I DDELL WITH MUT5 FROM YOS 

H 

122 

C 


H 

123 

16 

IF (TT5.LF. 1 500 . ) GO TO 17 

H 

124 


CALL DISCOT (TT5.PT5.TABTY.T ABNUY « T APY .11 *52 .4 »MUT5 ) 

H 

125 


GO TO 18 

H 

126 

1 7 

MUT5= 1 . 4 62 E -6* SORT ( TT5 )/( 1 .+112./TT5) 

H 

127 

1 B 

QT5FR=.760* ( (RHOW*MUw )**• 1 ) * ( ( RHOT 5*MUT5 ) ** * 4 )*(HT5-HW )*SQRT ( BO ) / ( 

_H_ 

128 


1 (PR**0.6)*2.**(SPH/2. 1 ) 

H 

129 


IF (ITEST.E0.4) GO To 20 

H 

130 

r 


H 

131 

c 

TOLQT REPRESENTS TOLERANCE OF ITERATION ON QT5 

H~ 

132 

c 


H 

133 


IF ( ABS ( 1 .-QT5M/QT5FR ) .LE. TOLQT > GO TO 20 

H 

134 


IF ( mM • EO • 1 5 ) GO TO 19 

H 

135 


MM=MM+ 1 

H 

1 36 


HT £ -= (HT5-HW )*QT5M/OTBFR + H\a/ 

H 

137 


IF (cAV.FO* 1 • ) GO TO 25 

_H_ 

138 


CALL SEARCH (PT5.RH0T5.HT5. ST5R .TT5.AT5.ZT5* GAMET5 . ZSTART5 » I SP* 1 ) 

H_ 

1 39 


GO TO 16 

H 

140 

^5 

CALL SAVF (PT5.RH0T5.HT5.ST5R.TT5.AT5.ZT5.GAMET5.3) 

H 

141 


GO To 1 6 

H 

142 

1 9 

PRINT 21 

H 

143 


I COOP* = 1 

H 

1 44 
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I 


20 

RETURN 


H 145 

c 



H 146 

21 

FORMAT ( 5 1 H ITERATIONS ON HT IN SECTION 3 OF QDOT EXCEED 

LIMIT) 

H 147 


END 


H 148 


SUBROUT INF SAVE (Pi R hO , H ♦ SR t T ♦ AM , Z , GAME ♦ K ) 


I 

1 

C 



I 

2 

c 

SAVE OBTAINS THERMODYNAMIC PROPERTIES FOR REAL AIR 


I 

3 

c 

IS BASED ON CURVE FIT EXPRESSIONS OF AEDC-TDR-63- 1 38 


I 

4 

c 

EXPRESSIONS OF AEOC-TDR-63- 1 38 APPLICABLE FOR T = 90 TO 

15000 

I 

5 

c 



I 

6 

c 

MAXIMUM PERCENT ERRORS- T=2000 TO 15000. AND P=lE+4 TO 

' 1E+6 

I 

7 

c 



I 

8 

c. 

RHO H T A 2 

game 

I 

9 

r 

2.42 1.96 2.24 2.78 0.75 

5.68 

I 

10 

r 



I 

1 1 

c 

INPUTS ARE PRESSURE fN/SQ METER) AND- 


I 

12 

r 

(1 ) ENTROPY, Q/R (K=l ) 


I 

13 

C 

(2) DENSITY. KG/CUBIC METER (K = 2) 


I 

14 

c 

(3) ENTHALPY, SQ METER/SQ SEC (K=3) 


I 

15 

c 



I 

16 

c 

ALSO, INPUTS DENSITY AND ENTHALPY ARE INCLUDED <K = 4 ) 


I 

1 7 

c 



I 

18 


DIMENSION T ABSR ( 6 ) , TABR(6)* TABHC6) 


I 

19 


DIMENSION T ABPM (13). T ABHM (13). T ABSRM (13), PM ( 1 3 ) 


I 

20 


NN=0 


I 

21 


MM = 0 


I 

22 


IF (K.NE.4) GO TO 3 


I 

23 


CON5=*03 


I 

24 


no 2 J=1 , 1 3 


I 

25 


PM ( J ) =RHO*H*CON5 


I 

26 


TABPm( J)=PM( J) 


I 

27 


P = PM ( J ) 


I 

28 


GO TO 3 


I 

29 

1 

TABHM ( J ) =HA 


I 

30 


T ABSRM ( J)=SR 


I 

31 


C0N5=C0N5+. 03 


I 

32 


MM = 0 


I 

33 


NN = 0 


I 

34 

2 

CONTINUE 


I 

35 


CALL FTLUP ( H , P , 2 . 1 3 , T ABHM , T ABPM ) 


I 

36 


CALL FTLUP ( H , SR , 2 , 1 3 , TABHM , TABSRM ) 


1 

37 


MM = 3 


I 

38 

3 

PL0G=AL0G1 0 (P/1 .0132^E+5) 


I 

39 


A=PLOG*PL0G 


I 

40 
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CJ1 



C=A*PLOG 

I 

-4-1- 


IF (K-.FQ. 1 ) GO TO 5 

I 

42 


IF (K.EQ.4.AND.MM.EQ.3) GO TO 5 

I 

43 


SRUP= 1 42 • 

1 

44 


SRLOw= 14* 

I 

45 


SR = ( 5RUP-SRL0W )/2. + l 4* 

I 

46 

4 

DFL8r= (SRUP-SRLOW )/2. 

I 

47 


IF ( NN • EQ • 0 ) GO TO 5 

I 

48 


SR-SQUP— DFL8R 

I 

49 

5 

SRLOG= ALOG 1 0 ( SR ) 

I 

50 


R=SRlOG*SRLOG 

I 

51 


D=R*oRLOG 

I 

52 


X 1 5 = — 39 • 1 442+83 •0558*SRL0G- 38. 204 2*SRL0G*SRL0G 

I 

53 


XI 51 =-l 0.*< PLOG-X1 5 ) 

I 

54 


IF (X151-40*) 7,6*6 

I 

55 

6 

T 1 5 = 0*0 

I 

56 


GO TO 1 0 

I 

57 ... 

7 

IF (X151+40.) 8,9,9 

1 

58 

8 

T 1 5= 1 . 0 

I 

59 


GO TO 10 

I 

60 

9 

T 1 5= 1 • / ( 1 • +FXP ( X 1 5 1 ) ) 

I 

61 

1 * 

IF (x.EQ.3. AND. MM. NF. 2 ) GO TO 36 

I 

62 


IF (K-.EQ.2. AND. MM. EQ. 2 ) GO TO 36 

I 

63 


IF (k.FQ.4* AND.MM.FQ.3 ) GO TO 58 

I 

64 


IF ('<'. cr Q.4.AND.MM,FO,2 ) GO TO 36 

I 

65 

r 


I 

66 

c 

COMPUTING RHO AS A FUNCTION OF P AND S/R 

I 

67 

r 


I 

68 

1 1 

XR1 2=- 1 6*5527+57.45*gRLOG-30 .8036*8 

I 

69 


XR23=499 • 544—938 ,91 *e,RLOG+609» 028*B— 1 35.995*D 

I 

70 


XR34=360*50 7— 634 *538*SR!_0G + 389 • 1 74*B-82 • 4653*D 

I 

71 


XR45=489. 628-458, 5*SRL0G+ 1 06 » 25*B 

I 

72 


XRl 21 =- 1 0.* (PLOG-XR1 9 ) 

I 

73 


XR231 = - 1 O.tt (PLOG-^<R2-3 ) 

I 

74 


XR34 1 = -10.* (p L aG-XR34 ) 

I 

75 


XR45 1 =-10.* (PL0G-XR4R ) 

I 

76 


IF ( y R 1 21 —40 • ) 12*15,15 

r 

77 

1 2 

IF (xRl 21 +40. ) 13,14,14 

i 

78 

1 8 

TR 12=1 • 0 

i 

79 


GO To 16 

i 

80 

1 4 

T R 1 2 = 1 . / ( 1 ,+EXP (XRl 2] ) ) 

i 

81 


GO TO 1 6 

i 

82 

i_5 

TR 1 2 = 0 . 0 

i 

83 

16 

IF (XR231-40*) 17,20,80 

j_ 

84 
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1 


1 7 

IF (XR231+4C.) 18,19,19 

I 85 

1 R 

TR23=1 .0 

I 86 


GO TO p\ 

I 87 

1 9 

TR23 = 1 . / ( 1 . +EXP (XR231 ) ) 

I 88 


GO To 21 

I 89 

20 

T9?3 = 9 . 0 

I 90 

?! 

IF (XR341-40.) 22,25,25 

I 91 

22 

IF (XR341+40.) 23,24,24 

I 92 

23 

TR34 = 1 • 

I 93 


GO TO 26 

I 94 

24 

TR34=1 ./(l . + EXP ( XR34 1 ) ) 

I 95 


GO TO 26 

I 96 


TR34=0 * 0 

I 97 

26 

IF (XR451-40.) 27,30,30 

I 98 

27 

IF (XR451+40.) 28,29,29 

I 99 

28 

TR45=1 .0 

I 100 


GO TO 31 

I 101 

29 

TR45=1 ./( 1 . +EXP(XR45] ) ) 

I 102 


GO TO 31 

I 103 

30 

TR45=0.0 

I 1 04 

31 

RHCLl = 1 5.951 867-0. 00?28295*PLOG- 15. 994242*SRL0G+. 0065 187267* A+. 530 

I 105 


1 7968p*PL0G*SRL0G+3. 1 75974*6 

I 106 


RHCL?= 1 541.1 666-63 • 9 7035*PI_OG-29 93 . 1 662 * SRLOG+ , 935 437*A+84 . 30375*5 

I 107 


1 RL0G*PL0G+1 938.7061 *p- . 0047460 1 6*C- . 6 1 284 04*A*SRLOG-27 . 422666*B*PL 

I 1 08 


20G-M °.08R1 *D 

I 109 


RHCL3=427.4745— 18 . 1 2 A622*PL0G-765 . 4 7626*SRL0G+ . 29343 1 69* A+22 . 92687 

I 110 


1 7*P(_OG*SR|_OG+4 56.71 7*B- . 00 1 7033404*C- . 1 8 0683 09* A *SRLOG-6 • 9 1 4 36 1 7*B 

I 111 


2* D L0G-9 I .131851*0 

I 112 


RHCL4 = 206. 23144-8 .2270278*PLOG-329 • 5465*SRL0G+ . 1 32 4 1 9 1 * A + 9 . 8884 1 65 

I 113 


1 *PL0G*SRL0G+1 75*03931 *B-.001 01 784 54*C- . 0765437 1 *A*SRL0G-2 . 6920 1 44* 

I 114 


26-K-PLOG-3] .237834*0 

I 115 


RHCL6=- 399. 52358+ 1 2 . R<^9477*PL0G+4 11.641 44*SRL0G- . 0 976949 1 9*A-6. 220 

I 116 


1 4477*PLOG*RRLOG-1 06.6733*R 

I 117 


RHCAL=RHCL1 + (RHCL2-RHCL1 )*TR12+ (RHCL3-RHCL2 )*TR23+ ( RHCL4-RHCL3 ) *TR 

I 1 18 


1 34+ ( 9HCL5-RHCL4 )*tR4k 

I 119 


RH 1 5=- 79. 282533+6 .35 770 78*PL0G+1 79. 2272 1 *SRL0G- . 1 2607Q98*A-8 . 40 1 3 1 

I 120 


1 22*P|_0G*SRL0G-1 29 .95?69*B+. 001 00 37437*0 + . 094 1 855 1 1 *A*SRL0G+3. 12569 

I 121 


266*PL0G*B+ 30.203862*9 

I 122 


RHCAl = RH1 5+ ( RHCAL-RHi 5 ) *T 1 5 

I 123 


RHO A = ( 1 0.**RH0AL )*1 .79233 

I 124 


IF (k.FO. 1 ) GO TO 36 

I 125 


I F (V.e0«3) GO TO 58 

I 126 


IF ( k'.EQ • 2 « AMD « MM. EQ. 1 ) GO TO 56 

I 127 


IF (^»EQi4. AND.MM.FQ, 1 ) GO TO 56 


Oi 

CJi 


I 128 
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oi 
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c 

CONVERGENCE TEST FOR K = 2 

I 

130 

r. 


I 

131 


IF ( ABS ( 1 .-RHQ/RHOA ) .LE. .001 ) GO TO 36 

I 

132 


NN=NM+1 

I 

133 


IF (PHO.GT.RHOA ) GO TO 34 

I 

134 

32 

SRLO\./ = SR 

I 

1 35 


SRUP=SR+DFLSR • 

1 

136 


IF ( nELSR.GT » 1 • ) GO TO 4 

I 

137 


TABSR ( 1 ) =SRLOW 

I 

138 


TABSR ( 6 ) =SRUP 

I 

139 


IF (K.EQ.2.0R.K.EQ.4 ) GO TO 33 

I 

140 


tr 

! 

1 

141 


GO TO 49 

I 

142 


N = ~? 

I 

143 


GO TO 64 

I 

1 44 

34 

SRUP=SR 

I 

145 


SRL0'»i = SR-DELSR 

I 

146 


1^ (OELSR.GT. 1 . ) GO TO 4 

I 

147 


TABSR ( 1 )=SRUP 

I 

148 


TABSP ( 6 ) =SRLOW 

I 

149 


IF (K.EQ.2.0R.K.EQ.4 ) GO TO 35 

I 

150 


1 

1 

1 

OJ 

li 

Z 

I 

151 


GO TO 49 

I 

152 

3S 

N = 2 

I 

153 


GO TO 64 

I 

154 

r 


I 

155 

" c 

COMPUTING ENTHALPY AS A FUNCTION OF P AND S/R 

I 

156 

r 


I 

157 

26 

IF (c;RLOG-l . 6 ) 37,37,38 

I 

158 

37 

HRCA|_= 1 2 .693869+5.397531 2*PL0G-48. 72 92 1 7*SRL0G-. 14 96152 1 *A-5. 87887 

I 

1 59 


1 74*Pl0G*SRL0G+48. 1 9278*B+ * 00090 1 44 1 32*C+. 091 151 473* A*SRL0G+ 1 .62828 

I 

_160_ 


2 29*P|_0G*B- 1 3.065267*0 

I 

161 


GO TO 48 

I 

162 

38 

IF ( GRL OG— 1 .76) 39 ,39, 46 

I 

163 

39 

HR22=- 1 56.37194+6 • 6959228*PL0G + 269 • 93097*SRL0G- .0971 79965* A -7 .5379 

I 

164 


1 71 4*PL0G*SRL0G-1 52. 1 3866*B+ . 00057029937*C + . 058364795*A*SRLOG + 2 . 1 59 

I 

165 


2 2 755*PL0G*B+28. 940926*0 

I 

166 


HR21 =-84 .008522+2 .5761 31 8*PL0G+ 1 07 . 06 1 98* SRLOG- . 0 1 4352904*A-1 .5313 

I_ 

167 


1 i 94 *dlOG*2RLOG-32 .31 6439*R 

I 

168 


XH =- 61 . 2053+114.1 03*fRLOG-47.5532*B 

I 

169 _ 


XH1 =-10.* (PLOG-XH ) 

_JL 

170 


IF (VH1-40,) 40,43,4^ 

I 

171 

40 

IF ( V H 1 +40, ) 41 ,42,4? 

i 

172 


t 
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41 

TH= ! . 

I 

173 


GO TO 44 

I 

174 

42 

TH= 1 ./( 1 .+EXPIXH1 ) ) 

I 

175 


GO TO 44 

I 

176 

43 

TH=0 • 0 

I 

177 

44 

HPCAL=HR21 + (HR22-HR21 )*TH 

I 

178 


GO TO 48 

I 

179 

4 *3 

IF ( SRLOG— I • 92 ) 46.46.47 

I 

100 

46 

HRC Al = — 35 • 1 60671 + • 5366924*PI_0G+56 • 99585*SRl_0G- • 02266 1 358* A- • 484703 

I 

181 


1 05*SRI_OG*PLOG-27.641 087*B+ . 00058568839*C+ . 0 1 6299962*A*SRL0G+ . 1 4073 

I 

182 


2606*R*PLOG+4. 712261 *D 

I 

183 


GO TO 48 

I 

184 

47 

HRCAl=- 1 1 4 • 94796+4. 004583*PLOG+1 80 • 08427*SRL0G- . 04 1 327787* A -4 • 0366 

I 

185 


1 535*PLOG*SRLOG-90 .76006*B+ . 00040320694*C+ . 024360248*A*SRLOG+ 1 .0462 

I 

186 


2299*PL0G*B+1 5.467804*0 

I 

187 

48 

HR 1 5 = 28 . 1 60664—2 . 2339873*PL0G— 59 . 053694*SRLOG+ . 054 97354 4* A+3 ,71 832 

I 

188 


15t*plOG*SRLOG+ 40.986503*B-.0004292698*C-. 040726332* A*SRLOG-l .37045 

I 

189 


205*PLOG*B-8. 253645*0 

I 

190 


HRCAl_ = HR 1 5+ (HRCAL-HR1 5 )*T1 5 

I 

191 


HA = ( 1 0.**HRCAL >*287.0 388 

I 

192 


IF (K.EQ.l .OP.K.EQ.2 ) GO TO 58 

I 

193 


IF ( K • EQ • 3 • AND • MM • EQ # 1 ) GO TO 52 

I 

194 


IF ( .< • E Q • 4 ) GO TO 1 

I 

195 

C 


I 

196 

r 

CONVERGENCE TEST FOR K=3 

I 

197 

C 


I 

198 


IF (ABS ( 1 .-H/HA > .LE. .001 ) GO TO 11 

I 

199 


NN=NN+1 

I 

200 


IF (HA.GT.H) GO TO 34 

I 

201 


GO TO 32 

I 

202 

r 


I 

203 

c 

interpolation for delsr less than i 

I 

204 

c 


I 

205 

49 

T ABH ( 1 )=HA 

I 

206 

50 

DELSR= ( TABSR (6 ) -TABSR ( 1 > >/5. 

I 

207 


DO 51 1=2.5 

I 

208 


T ABSR ( I )=TABSR( 1-1 ) +OELSR 

I 

209 

51 

CONTINUE 

I 

210 


IF ( • E Q « 2 • OR • K • f- Q • 4 ) GO TO 55 

I 

21 1 


DO 57 1=2.6 

I 

212 


MM = 1 

I 

213 


SR=TABSR( I ) 

I 

214 


GO TO 5 

I 

215 

^2 

T ABH ( I ) =HA 

I 

216 


CJ1 
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53 

CONTINUE 

J 

21 7 


CALL FTLUP <H,SR.N,6,TABH,TABSR ) 

I 

218 



I 

219 


GO TO 5 

I 

220 

54 

TABR ( 1 ) =RHO A 

! 

221 


GO TO 50 

I 

222 

55 

no 57 1=3,6 

I 

223 


MM = 1 

I 

224 


SR=TABSR< I ) 

I 

225 


GO TO 5 

I 

226 

^6 

TARR ( \ ) =RHOA 

1 

227 

57 

CONTINUE 

I 

228 


CALL FTLUP ( RHO , SR . N , 6 ♦ TABR » TABSR ) 

I 

229 


MM = 2 

1 

230 


GO TO 5 

I 

231 

C 

- 

1 

232 

C 

COMPUTING Z 

I 

233 

c 


I 

234 

5Q 

XZ12=62. 91-41 .5*SRL0G 

1 

235 


XZ23=72. 945-45. 75*SRlOG : 

I 

236 


XZ 34 = 65 • 75— .77 • 5*SRLO0 

I 

237 


X745-62 .92-32 .0*SRLOG 

I 

238 


X71 21 =- 1 o .* (PL0G-X71 P ) 

I 

239 


X2231 =-10.* (PL0G-XZ23 ) 

1 

240 


XZ341 =-10.* (PL0G-X734 ) 

I 

24 1 


XZ45! =-10.* (PL0G-XZ4R ) 

I 

242 


ZCAL?=51 9. 80374 3. 7535 14*PL0G-983. 90729* SRLOG+. 37296957* A+30. 0843 

I 

243 


1 79*PL0G*SRL0G+620 .04 1 68*B- . 002 1 648826*C- . 237 1 0079* A*SRL0G-9 . 496903 

I 

244 


2*PLOG*B- 1 29 . 7892 1 *D 

I 

245 


ZCAL3= 366 .40674—15.51 7444*PL0G— 647 , 42436*SRL0G+ . 1 870 1 758* A+ 1 8 . 0403 

I 

246 


1 83*PL0G*SRL0G+379 .59834*B-. 00087958438*C~. 1 0580 1 29*A*SRL0G-5 . 18882 

I 

247 


254*PL0G*B-73.504?6R*n 

1 

248 


ZCAL4 =51 6.07331 -16. 5o277*PL0G-808.49823*SRL0G+ .071 256235* A+l 6.5268 

I 

249 


1 1 3*PL0G*SRL0G+41 8.4534 1*B+. 00094 1 83347*C-. 01 972781 7*A*SRL0G-3.9 948 

I 

250 


2a06-M-PL0G*R-7i .038921 *o 

T 

251 


IF (XZ121-40.) 59,62,62 

I 

252 

59 

IF ( yZ 1 2 1 +40 • ) 60,61 ,61 

I 

253 

62 

TZ12= 1 • 

1 

254 


GO TO 63 

r 

255 

61 

TZ12=1 •/( 1 • +EXP ( XZ 1 2 1 ) ) 

i 

256 


GO TO 63 

i 

257 

63 

TZ12=0.0 

i 

258 

63 

IF (XZ231-40.) 64,67,67 

l 

259 

64 

TF (XZ231+40.) 65,66,66 

i 

260 
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6^ 

T Z 2 3 = 1 . 

I 

261 


GO TO 68 

I 

262 

66 

TZ23=1 ./Cl .+FXPCXZ231 ) > 

I 

263 


GO TO 68 

1 

264 

67 

T723=n.O 

I 

265 

68 

IF ( XZ34 1 —40 • > 69,72,72 

I 

266 

69 

IF (XZ341+40. ) 70,71 ,71 

I 

267 

70 

TZ 34 = 1 . 

I 

268 


GO TO 73 

I 

269 

7 i 

TZ34= 1 • / < 1 .+EXP (XZ34 1 )) 

I 

270 


GO TO 73 

I 

271 

7 P 

TZ34=0.0 

I 

272 

77 

IF (XZ451-40. ) 74,77,77 

I 

273 

74 

IF ( XZ4 5 1 +4 0 « ) 78,76,76 

I 

274 

77 

TZ4F=1 . 

I 

275 


GO To 78 

I 

276 

76 

TZ45= 1 ./(l • +FXP ( XZ45 1 ) ) 

I 

277 


GO To 78 

I 

278 

77 

TZ45=0 • 0 

I 

279 

78 

ZCAL=1 .0+ (ZCAL2-1 . ) *TZ 1 2 + ( ZC AL3-ZC AL2 1 *TZ23 + ( ZCAL4-ZCAL3 ) #TZ34+ (4. 

I 

280 


1 0-ZCAL4 )*TZ4fJ 

I 

281 


ZCA|_=1 . + (ZCAL-1 . >*T1F 

I 

282 


7 = ZC A L 

I 

283 

r 


I 

284 

r 

COMPUTING T ( DEG K) 

I 

285 

c 


I 

286 


W0 = 28 # 967 

I 

287 


RUNl V=831 4. 34 

I 

288 


IF CK-.EQ.2.0R.K.EQ.4 ) GO TO 79 

I 

289 


RHO=RHO A 

I 

290 

79 

T=P*wO/ (RH0*RUN1V*Z ) 

I 

291 

c 


I 

292 

c 

COMPUTING A (M/SEC) 

I 

293 

c 


1 

294 


IF (t-2100. ) 80 <80(81 

I 

295 


IF (T-I500. ) 82 , 82 , 8 1 

I 

296 

81 

IF (PLOG+1 . ) 83.83,8? 

I 

297 

8? 

CONI = SORT ( T/273. 1 5 ) 

I 

298 


AOAO = - « 0753 8 08+ CON 1 * t 1 . 1 2644- . 0552696*CON 1 ) 

I 

299 


AM=371 .31 1 5*A0A0 

I 

300 


GO TO 104 

I 

301 

83 

XA 12=635.054-1 220 .46*SR|_0G+8 03. 882*B-1 80 .84 5*D 

I 

302 


X A23=373. 702-663. 358*SRL0G+4 08. 854 *B-86. 8 0B6*D 

I 

303 


XA 34 = 1 703.78-2602 .97*SRL0G+1 337. 93*B-231 .422*D 

I 

304 


CJ1 
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XA 22= 1 043 ,37- 1 820 .34*SRL0G+1 076 • 36*B-2 1 5 • 44 5*D 

I 

305 


XA12t =-10.* (PLPG-XA 1 2 ) 

I 

306 


XA231 =- 1 0 .* (PLOG-XA2F ) 

I 

307 


XA341 = - 1 0 .* ( PL0G-XA34 ) 

I 

308 


XA22 1 = -l 0.* (PL0G-XA22 ) 

I 

309 


A 1 =-4409.6241+1 96.82 ?59*PL0G+8746. 4 634*SRL0G-3. 1650299* A-262. 32947 

I 

310 


1 *PL0G* SRLOG- 5786. 449*B+. 02 00 04 1 86*C+2 • 1 4 2982 5* A*SRL0G+87 • 58 9029 *PL 

I 

31 1 


20G*B+1277. 6718*0 

I 

312 


A2 1 =— 1 8 1 4 • 5 1 1 7 + 86 .0 960 78*PL0G+33 1 5 • 6099*SRl_0G- 1 • 7593034* A— 1 07.2534 

I 

313 


1 *PLOG*SRLOG-2023. 20 1 *B+ . 0 1 62 87679*C+ 1 . 1 398 1 34* A*SRL0G+33 . 659607*PL 

I 

314 


20G *B + 4 1 3.41945*0 

I 

315 


A22=?65 1 .2944-81 • 40 C; 596*PLOG-3099 • 0064*SRLOG+ • 69752668* A+ 48 • 062596 

I 

316 


1 *PL0G*SRLGG+907.7088o*B 

I 

317 


IF (XA221-40.) 85,84,84 

I 

318 

8 4 

TA?2=0. 0 

I 

319 


GO T n 88 

I 

320 

8F 

IF (XA221+40.) 86,86,87 

I 

321 

88 

TA2P=1 .0 

I 

322 


GO TO 88 

I 

323 

87 

TA22=l./(l.+EXP(XA22l ) ) 

I 

324 

88 

A2=A? 1 + ( A22-A21 >*TA22 

I 

325 


A 3 = -321 7.80 37+1 95 . 34964*PL0G + 5348 . 2 1 43*SRL0G-4 . 6268475* A- 22 1 . 1 2705 

I 

326 


1*PLOG*SRLOG-2970.8649*B+. 0446 14358*0+2.70791 77*A*SRL0G+63 . 042803*P 

I 

327 


2L0G*R+553 . 1 2007*0 

I 

328 


44=16976.939-476. 1 0242*PLOG- 1 7445 . 3 1 5*SRL0G+3 . 6534 057* A +246 . 4 1 1 25* 

I 

329 


IPLOG-h-SRLOG+4486.31 1 8*B 

I 

330 


IF ( X A 1 2 1 —40 . ) 90,89,89 

I 

331 

«Q 

TAlPrO.O 

I 

332 


GO Tn 93 

I 

333 

QO 

IF (yA 1 21 +40. ) 91,9 1,92 

I 

334 

51 

T A 1 2= 1 .0 

I 

335 


GO TO 93 

I 

336 

92 

T A 1 2 = 1 . / ( 1 . +FXP ( X A l 2 1 1 ) 

I 

337 

93 

I ^ (XA231-40.) 95,94,94 

I 

338 

94 

TA2?=0. 0 

I 

339 


GO TO 98 

I 

340 

cTcj 

IF (XA231+40.) 96.96,97 

I 

34 1 


T A 2 3 = 1 . 0 

I 

342 


GO TO 98 

I 

343 

97~ ” 

” TA'^'=17/TlT+FXP <x'a?31 ) ) " ~ ~ ~~ 

f 

344 

““"98 

IF" ( X A34"T— 4*6V ) 100 ,99 ,99 

I 

345 

Qa 

TA34=0.0 

I 

346 


GO TO 1 03 

I 

347 

fn 6“ ' 

IF ( X A 34 1 +4 0 . ) 101, ini, 102 

I 

348 
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101 

T A 34= 1 • 0 


I 

349 


GO TO 103 


I 

350 

102 

TA34=1 •/( 1 • + EXP ( X A34 1 ) ) 


I 

351 

1 03 

AOAO= A 1 + t A2-A1 ) *TA 1 2 + ( A3-A2 ) *TA23+ ( A4-A3 ) *TA34 


I 

352 


AM=331 .31 15*A0A0 


I 

353 

C 



I 

354 

C 

COMPUTING GAMF 


I 

355 

r 



I 

356 

1 04 

GAME=W0*AM**2/ (RlJNl V*Z*T ) 


I 

357 


IF (k.EG. 2) GO TO lOq 


I 

358 


IF (K.EQ.3.0R.K.EQ.4 ) GO TO 106 


I 

359 


< 

I 

II 
I 


I 

360 


GO TO 106 


I 

361 

1 05 

H = HA 


I 

362 

106 

RFTURN 


I 

363 


END 


I 

364- 


SUBROUT INE SEARCH (P ♦RHOtHi ♦ SOR i T1 • A 1 iZl « GAM * ZS « ISP *K ) 

J 

1 

c 



J 

2 

c 

K=1 CORRESPONDS TO INPUTS P AND HI 


J 

3 

c 

K= 2 CORRESPONDS TO INPUTS P AND RHO 


J 

4 

c 

K=3 CORRESPONDS TO I mPUTS RHO AND HI 


J 

5 

r 



J 

6 


DIMENSION G(4)* Yl(4)* Y2(4), Y3(4)i Y4(4), Y5(4) 

t Y6 ( 4 ) . Y7 ( 4 ) 

J 

7 


DIMENSION I COUNT (25). JFLAG(25). Y ( 9 . 1 50 ) . P(25). 

RHO ( 25 ) 

J 

8 


DIMENSION SAVEH ( 25 » 4 ) » SAVER<25.4), SAVET(25.4). 

SAVE A ( 25 * 4 ) 

J 

9 


DIMENSION S AVEZ (25.4). SAVES(25*4). SAVEG(25*4)« 

SAVEZS (25.4 ) 

J 

10 


DIMENSION S AVEP ( 25 . 4 ) 


J 

1 1 


DIMENSION HI (25 ) t Tl(25). Al<25). Zl(25). SOR ( 25 ) 

« GAM (25). ZS ( 25 ) 

J 

12 


DIMENSION T ABT ( 1 50 ) < TABR<150>. TABP ( 1 50 ) . TABHU50) 

J 

13 


DIMENSION T ABA ( 1 50 ) . TABZU50). TABG < 1 50 ) . TABZSU50) 

J 

14 


DO 1 1=1 • ISP 


J 

15 


I COUNT ( I ) = 1 


J 

16 


JFLAO ( I) =0 


J 

17 

1 

CONTINUE 


J 

18 


JUMP=0 


J 

19 


IT=8 


J 

20 


rewind it 


J 

21 

2 

READ (IT) X,NV« ( (Y< I ,L ) * 1 = 1 *9) »L=1 ,NV) 


J 

22 


IF (PNDFILE IT) 3«9 


J 

23 

3 

CONTINUE 


J 

24 


WRITF (6.40) 


J 

25 


DO 8 I =1 . ISP 


J 

26 


IF ( JFL AG (I ) • EQ .0 ) GO TO 4 


J 

27 


GO TO 8 


J 

28 


03 
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4 

CONTINUE 




J 

29 


IF (kr.EQ. 1 ) GO TO 5 




J 

30 


IF (x.EQ.2) GO TO 6 




J 

31 


P< I 1=0. 




J 

32 


GO TO 7 




J 

33 

5 

RHO ( I ) =0 • 




J 

34 


GO TO 7 




J 

35 

6 

HI (I 1 = 0. 




J 

36 

7 

T 1 (11=0. 




J 

37 


ai ( i >=o. 




J 

38 


71 ( I 1=0. 




J 

39 


SOR( 11=0. 




J 

40 


GAM ( n =0. 




J 

41 


zs c n=o. 




J 

42 

B 

CONTINUF 




J 

43 


GO TO 09 




J 

44 

o 

CONTINUE 




J 

45 


O0 3R J=1 , ISP 




J 

46 


IF ( JFLAG ( J ) .EQt M GO TO 38 




J 

47 


NN= I COUNT ( J ) 




J 

48 


IP (K.EQ. 1 ) GO TO 10 




J 

49 


IF (vr.FO.7 ) GO TO 11 




J 

50 


HH = A(_OGlO(Hl ( J)/287.n24F) 




J 

51 


RR = AL0G1 0 (RHO ( J 1/1 .20 1 489 ) 




J 

52 


GO TO 1 3 




J 

53 

1 0 

PP = Ai_OG10 (P( J )/l .01 3PFE+5 ) 




J 

54 


HH=AL0G10(H1 (J)/?87.o245) 




J 

55 


GO T P IP 




J 

56 

1 1 

PP = A|_OG 1 0 (P ( J 1/1 . 01 325E+5 ) 




J 

57 


RR = A|_OG 1 0 (RHO ( J )/l .291489) 




J 

56 

12 

IF ( ( PP — Y ( 3 i 1 ) )* (PP-y (3.NV ) ) .LT.O. 

) GO 

T0_ 

15 

J 

59 


GO TO 14 




J 

60 

1 3 

I F ( (RR-Y (? , 1 ) 1* (RR-y (2.NV ) ) .LT.O. 

) GO 

"to - 

1 5 

J 

61 

t 4 

SAVFD ( J, 1 ) =0. 




J 

62 


S A VFh ( J n 1 ) =0 * 




J 

63 


S AV^R ( Ji 1 1=0. 




J 

64 


SM/FT ( J, 1 ) = 0 • 




J 

65 


SAVFA(J« 1 ) = 0 • 




J 

66 


SAVF7 ( J . 1 ) =0. 




J 

67 


SAVE*; ( J , 1 )=0. 




J 

68 


SAVFG ( J* 1 ) =0. 




J 

69 


SAVEZS ( J* 1 ) =0. 




J 

70 


NN = 2 




J 

71 


GO TO p 4 j 72 
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1 s 

DO 16 T =1 t MV 

J 

73 


TABTf I)=Y(1, 

I ) 

J 

74 


TAPP ( 1 >=Y <?_i 

1 ) 

J 

75 


T ABP ( I )=V(.1, 

I ) 

J 

76 


TABH ( I I =Y (4t 

I > 

J 

77 


T AEG ( 1 > =Y ( 6 , 

I > 

J 

78 


TABA (I ) =Y (6 t 

I ) 

J 

79 


T ABZ ( I )=Y(7, 

I ) 

J 

80 


T ABZs C I ) = Y ( 9 

t I ) 

J 

81 

16 

CONTINUE 


J 

82 


IF (kr.EQ.3) 

GO TO 1 7 

J 

83 


CALL D I SCOT 

(PP tPP, TABP, TABH tTABH t-1 30tNV tO t ANSI ) 

J 

84 


CALL DISCOT 

( pp ,pp, TABP, TARR t TABR,-1 30 tNV t 0 t ANS2 ) 

J 

85 


CALL D I SCOT 

(PP ,PP, TAPP, TART . TABT ,-l 30, NV *0 , ANS3 ) 

J 

86 


CALL DISCOT 

(PP.PP, TABP. TABA .TABA , - 1 30 , NV , 0 . ANS4 ) 

J 

87 


CALL D I SCOT 

(PP,PP,TABP,TABZ.TABZ.-130.NV.0. ANS5) 

J 

88 


CALL D I SCOT 

(PPtPP, TABPt TA8G tTABG t-1 30tNVt 0 t ANS6 ) 

J 

89 


CALL DISCOT 

(PP ,pp, TABPt TARZStTABZS t-1 30 tNVt Ot ANS7 ) 

J 

90 


SAVFM JtNN)= 

X 

J 

91 


SAVEH ( NN ) = 

ansi 

J 

92 


SAVER ( J t NN ) = 

ANS2 

J 

93 


SAVFT ( J t NN ) = 

ANS3 

J 

94 


SAVEA ( J,NN> = 

ANSA 

J 

95 


SAVF7( JtNN)= 

ANS5 

J 

96 


savfgi j,nni= 

A NS 6 

J 

97 


SAVF7S ( Jf NN ) 

= ANS7 

J 

98 


CO TO 18 


J 

99 

1 7 

CALL DISCOT 

(RRiRRtTABRtTABH iTABHi-130»NV (OiANSl ) 

J 

100 


CALL DISCOT 

(RR .RR, TABR , TABT tTABT i-130tNVt0t ANS3 ) 

J 

101 


CALL D I SCOT 

(RR.RR, TAt-R, TABA .TABA , - 1 30 . NV , 0 , ANSA ) 

J 

102 


CALL DISCOT 

(RR.RR,TABR.TABZ,TABZ,-130.NV.0»ANS5) 

J 

103 


CALL D I SCOT 

(RR .RR, TABR . TABG .TABG.-l 30.NV .0. ANS6 ) 

J 

1 04 


CALL D I SCOT 

( RR, RR, T ABR, T ABZS, T ABZS. -1 30, NV.0.ANS7) 

J 

1 05 


CALL D I SCOT 

( RRt RR tT ABR t TABP t TABP t-1 30 tNVtQt ANS8 ) 

J 

106 


SAVFS( JtNN)= 

X 

J 

107 


SAVEHI J.NN)= 

ANSI 

J 

1 08 


SAVEt ( J'NN) a 

ANS3 

J 

1 09 


SAVEA ( JtNN ) = 

ANS4 

J 

no 


S A VE7 ( J t NN ) = 

ANS5 

J 

1 1 1 


SAVES (JtNN) = 

ANS6 

J 

1 12 


SAVE7S ( JtNN ) 

= ANS7 

J 

113 


SAVEP ( JtNN ) = 

ANS8 

J 

114 

1 8 

IF (<.FQ. n 

GO TO 20 

J 

1 15 


IF <kr.FQ*2> 

GO TO 21 

J 

116 
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IF ( SAVEH ( J * NN ) • GT • HH ) GO 

TO 

25 


J 

1 1 7 


IF ( \tN • EQ • 3 ) GO TO IQ 




J 

1 18 


NN=NN+ 1 




J 

119 


GO TO 24 




J 

1 20 

1 9 

SAVFP ( J t 1 ) = SAVFP ( J , 2 ) 




J 

121 


SAVEH ( J* 1 >=SAVFH< J,2 ) 




J 

122 


FA VET ( Jt 1 ) =S AVET ( J»2 ) 




J 

123 


SAVE A ( Ji 1 ) =S AVFA < J , 2 ) 




J 

124 


SAVE/ ( J t 1 ) =SAVEZ ( J » 2 1 




J 

125 


SAVFF( J , 1 ) = SAVES ( Ji2 ) 




J 

126 


3 A VEG ( J « 1 ) =SAVEG ( J*2 ) 




J 

127 


SA VE/S ( J i 1 ) =SAVEZS( J,2 > 




J 

128 


SAVFP ( J* 2 ) =SAVFP (J,3) 




J 

1 29 


SAVEH < J t 2 1 =SAVEH ( J , 3 ) 




J 

130 


SAVET ( J«2)=SAVFT( J,3 ) 




J 

131 


SAVFA ( J*2 ) = SAVE A ( J,3 ) 




J 

1 32 


SAVE/ ( Jt2 ) = S A VE Z ( J , 2 1 




J 

1 33 


SAVES ( J i 2 ) =SAVES ( J«3 ) 




J 

134 


SAVFO ( J,2 ) = SA VFG ( J , 3 i 




J 

135 


SAVE/S ( Ji 2 )=SAVEZS ( J » 3 ) 




J 

1 36 


GO TO 24 




J 

137 

20 

IF (SAVEH (Ji NN ) • GT • HH ) GO 

TO 

25 


J 

1 38 


GO T n 22 




j 

139 

21 

IF ( SAVER ( JtNN ) .LT.RR ) GO 

TO 

25 


J 

140 

22 

IF (MN.FQ.3) GO TO 2 7 




J 

141 


NN=NM+1 




J 

142 


GO TO 24 




J 

143 

22 

SAVER (Jtl ) =SAVER ( J » 2 ) 




J 

1 44 


SAVEH (Jil ) = S A V E H ( J , 2 ) 




J 

145 


SAVET <J»1 )=FAVET(J,2) 




J 

1 46 


SAVE A ( J « 1 ) = SAVE A ( J , 2 ) 




J 

147 


SA VE7 ( 1 )=SAVEZ(J,2) 




J 

148 


SAVES ( J t 1 ) =S AVES (Ji2) 




J 

149 


S A VEG ( J t 1 ) = SAVE' - ( Ji? ) 




J 

1 50 


SAVS7S ( J t 1 )=SAVEZS( J,2) 




J 

151 


SAVER (Jt2 ) = SAVER ( J t 3 ) 




J 

152 


SAVEH ( J«2 )= SAVEH ( J, 3 ) 




J 

1 53 


SAVFT ( Jn2 ) =5 A VET C J, 3 ) 




J_ 

1 54 


SAVFA ( Jt? ) =SAVFA ( J ,3 ) 




J 

155 


SAVE 7 ( J , 2 ) =SA VFZ ( J , 3 ) 




J 

1 56 


SAVES ( Jt2 ) =SAVES ( J, 3 ) 




J 

157 


SA VEG ( J t 2 ) = S A VEG ( J, 3 ) 




J 

158 


SAVE/S ( J .2 ) =SAVEZS ( J t 3 ) 




J 

159 

24 

I COUNT ( J)=NN 




J 

160 
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GO TO 38 

J 

161 


IF ( NjN • EO • 4 ) GO TO 26 

J 

162 


NN=NN+1 

J 

163 


I COUNT ( J)=NN 

J 

164 


GO To 38 

J 

165 

26 

JFL Ao ( J ) = l 

J 

166 


00 3 n M = 1 ,4 

J 

167 


IF (x-.FQ. 1 ) GO TO 27 

J 

168 


IF Oc.EQ.2) GO TO 28 

J 

169 


G ( M ) =SAVEH ( J » M ) 

J 

170 


Y 1 (M)=SAVEP<J*M> 

J 

171 


GO TO 29 

J 

172 

27 

G(M)=SAVEH(J*M) 

J 

173 


Y I CM ) =SAVER ( J .M ) 

J 

174 


GO TO 29 

J 

175 

28 

G ( M ) = S A VER ( J ♦ M ) 

J 

176 


Y 1 <M )=SAVFH( J,M ) 

J 

177 

29 

Y2 (M ) =S AVFT (J«M) 

J 

178 


Y3 CM ) =SAVEA (JiM) 

J 

179 


Y4 CM ) =SAVEZ I J(M > 

J 

180 


Y5 CM )=SAVES < J*M ) 

J 

181 


Y6 (M ) =SAVEG C J»M > 

J 

182 


Y7(M)=SAVEZS(J,M) 

J 

183 

30 

CONTINUE 

J 

184 


IF (kT.FQ. 1 ) GO TO 31 

J 

185 


IF UC.EQ.2) GO TO 33 

J 

186 


CALL INTRP (4,G»Y1 ,Hh.P) 

-J 

187 


GO TO 3? 

J 

188 

31 

CALL INTRP <4«G«Y1 ,HH«R> 

J 

189 

32 

CALL INTRP (4*G*Y2*HH«T) 

J 

190 


CALL INTRP (4*G.Y3,HhU A) 

J 

191 


CALL INTRP (4«G* Y4,HmZ) 

J 

192 


CALL INTRP ( 4 ♦ G « Y5 * HH « SR 1 ) 

J 

193 


CALL INTRP ( 4 » G » Y6 » HH * GAM 1 ) 

J 

194 


CALL INTRP (4,GfY7,HH*ZSl ) 

J 

195 


GO TO 34 

J 

196 

33 

CALL INTRP ( 4 * G * Y 1 * Rp * H ) 

J 

197 


CALL INTRP (4iGtY2i RR • T ) 

J 

198 


CALL INTRP (4iGtY3i RR t A ) 

J 

199 


CALL INTRP (4iG'V4i RR i Z ) 

J 

200 


CALL INTRP <4»G«Y5»RR«SR1 ) 

J 

201 


CALL INTRP ( 4 » G » Y6 « RP « GAM 1 ) 

J 

202 


CALL INTRP (4,G*Y7*RR«ZS1 ) 

J 

203 

34 

IF (k-.EQ.I ) GO TO 35 

J 

204 


a> 

U1 
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IF ( kT • EO • 2 ) 00 TO 36 

J 

205 


P ( J) = ( 1 0.**P )*1 .01 32=>E+5 

J 

206 


GO TO 37 

J 

207 

35 

RHO ( J ) = ( 1 0. **R ) *1 .291489 

J 

208 


GO TO 37 

J 

209 

36 

HI ( J ) = ( 1 0.**H )*287.0?4 5 

J 

210 

3 7 

T 1 ( J ) =T 

J 

21 1 


A 1 ( J ^ = A *33 1 *41 84 

J 

212 


71 ( J ) =Z 

J 

21 3 


SOR ( J ) = SR 1 

J 

214 


GAM ( J ) =GAM 1 

J 

215 


2S ( J )=ZS1 

J 

216 


JUMP= JUMP+1 

J 

217 


IF { JUMP. EO. ISP ) GO TO 39 

J 

218 

38 

CONTINUE 

J 

219 


GO TO 2 

J 

220 

39 

CONTINUE 

J 

221 


RFTUPN 

J 

222 

c 


J 

223 

40 

FORMAT ( 1 HI «60X. 7HWARN I NG//// ) 

J 

224 


FND 

J 

225- 


SUBROUT INF SLOW ( XX i7i IhJl » I T • N V * NERR i V * X ) 

K 

1 

C 

TAPE. IS WRITTEN WITH LINES OF CONSTANT XX 

K 

2 

c 

Z(I1) AND XX ARE INDfPENDENT VARIABLES 

K 

3 

c 

Z(J1) IS THE DEPENDENT variable 

K 

4 

c 

AK= +1. IE XX INCREASES MONOTON I CALL Y ON TAPE 

K 

5 

c 

AK= -1. IF XX DECREASES MONOTON I CALL Y ON TAPE 

K_ 

6 

c 

IT- TAPE UNIT 

K 

7 

c 

NV= NO. OF VARIABLES ON TAPE FOR EACH XX .(NOT GREATER THAN 9) 

K 

8 

c 

NO. OF POINTS FOR EACH XX NOT GREATER THAN 150 

K 

9 

c 

BEGIN EXECUTION 

K 

10 


DIMENSION X ( 4 ) * Y ( 4 , Q , 1 50 ) . Z<9), U(4), V(4), W(4>. NP(4) 


1 1 


COMMON I COUNT , I MET ( 2 ) . NP , AE A R . ME . ME 



K 

12 


RFAL MF.MF 

K 

13 


I COUNT = I COUNT + 1 

K 

14 


IE t J MET ( 1)) 3,1,3 

K 

15 

1 

BACKSPACE IT 

K 

1 6 


READ (IT) DUM 

K 

1 7 


rewind it 

K 

18 


DO 2 K = 1 .3 

K 

19 


READ (IT) X (K ) , J , ( ( Y (K , I ,L ) , I = 1 ,NV ) ,L= 1 , J ) 

< 

20 

2 

NP (K )~J 

K 

21 


XW = X (2 ) -X ( 1 ) 

K 

22 


AK= APS ( XW ) /XW 

K 

23 
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niRi =i . 

K 

24 


IMFTf 1 ) = 1 

K 

25 


xxx=xx 

K 

26 


NFRR=0 

K 

27 


I M = 3 

K 

28 


GO Tn 18 

K 

29 

3 

NEPR=0 

K 

30 

C 

EXCEPT FOR FIRST TIMr THROUGH 

K 

31 


IF < <XX-X(M1 ) )* <XX-X(M2 ) ) ) 25,25,4 

K 

32 

4 

TEMP= (XX-XXX )*AK 

K 

33 


DIR2=ABS(TEMP)/TFMP 

K 

34 


GO=nyRl *DIR2 

K 

35 


xxx=xx 

K 

36 


0 I R 1 =D I R2 

K 

37 


IF ( n I R2 ) 5.38,16 

K 

38 

r 


K 

39 

r 

negative direction 

K 

40 

5 

IF (GO) 6*35*7 

K 

41 

6 

BACKSPACE IT 

K 

42 


BACKSPACE IT 

K 

43 


BACKSPACE IT 

K 

44 


GO TO 9 

< 

45 

7 

j M= Im-1 

K 

46 


IF (TM) B* 8*9 

K 

47 

S 

ii 

£ 

K 

48 

9 

Ml =IM+1 

K 

49 


BACKSPACE IT 

K 

50 


BACKSPACF it 

K 

51 


IF (Ml -4 ) 1 1 1 1 1 * 1 0 

K 

52 

1 o 

Ml =1 

K 

53 

1 1 

M? = Mt +1 

K 

54 


IF ( m2— 4 ) 13.13,12 

K 

55 

12 

M2 = l 

K 

56 

13 

read (IT) X ( IM ) , J , ( < Y < IM, I ,L ) , I = 1 ,NV ) ,L= 1 , J ) 

K 

57 


NP(IM)=J 

K 

58 


IF ( fXX-X (Ml ) )* (XX-X (M2 ) ) ) 25,25,14 

K 

59 

1 4 

IF <X(M1 )-X(M2) ) 7,18,7 

K 

60 

C 

ERROR, VARIABLE OFF FRONT END OF TAPE 

K 

61 

1 =5 

CONTINUE 

< 

62 


NERR= 1 

K 

63 


GO TO 36 

K 

64 

r 


K 

65 

c 

POSITIVE DIRECTION 

K 

66 

1 6 

IF (GO) 17*35*18 

K 

67 


05 

-a 
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Oi 

CO 


1 7 

RFAn (IT) DUM 


K 

68 


RF AD (IT) HUM 


K 

69 


READ (IT) DUM 


K 

70 


GO TO 20 


K 

71 

18 

I M= IM+1 


K 

72 


IF (TM-4) 20,20*19 


K 

73 

1 9 

IM=1 


K 

74 

20 

Ml = I M- 1 


K 

75 


IF (Ml ) 21,21 ,22 


K 

76 

21 

M 1 =4 


K 

77 

2? 

M2=M1 -1 


K 

78 


IF (m2) 23,23,24 


K 

79 

23 

M2 = 4 


K 

80 

24 

READ (IT) X ( IM ) , J , ( ( Y ( I M , I ,L ) , I = 1 ,NV ) ,L= 1 , J ) 


K 

81 


NP ( IM ) = J 


K 

82 


IF < (XX-X (Ml > )* (XX-X (M2 ) > ) 25,25,18 


K 

83 

C. 



K 

84 

C 

TAPE SEARCH COMPLFTF , DO CROSS FOUR POINT 


K 

85 

36 

DO 34 K=1 ,4 


K 

86 


NPK=mP (K)- l 


K 

87 


DO 26 1=1 « NPK 


K 

88 


IF ( ( Y (K , I 1 , I )-Z ( I 1 ) ) * ( Y (K , I 1 , I + 1 )-Z ( I 1 ) ) ) 27,27,26 


K 

89 

26 

CONTINUE 


K 

90 


NFRR = 1 


K 

91 


GO TO 36 


K 

92 

27 

IF ( 1-1 ) 29,28,29 


K 

93 

28“ 

J = o 


K 

94 


GO To 33 


K 

95 

29 

IF { I —NPK ) 31,30,31 


K 

96 

3 ci 

J = NPk' — 3 


K 

97 


GO TO 32 


K 

98 

Vi 

J= 1-3 


K 

99 

32 

DO 33 L=1 ,4 


K 

io'b 


MX=L+J 


K 

101 


U ( L ) = Y ( K * I 1 *MX> 


K 

102 

33 

V ( L ) = Y ( K * J 1 ,MX> 


K 

1 03 


34 


IS 


CALL INTRP (4*U*V*Z ( ! 1 ) *W(K) ) 


CALL INTRP (4,X»W»XX,Z(J1 )) 


RFTUDN 


CONTINUE 


36 


nfrr= 1 


I F ( TMFT (2 ) ) 39* 37, 39 


3 T 


I MFT ( 2 ) = 1 


W°TT^ (6*40) XX* I 1 * 7 . ( I 1 ) ♦ J1 


K 1 04 


K 105 


K 106 


K 107 


K 1 08 


K 109 


K 110 


Kill 
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DO 3R I M= 1 ,4 

K 

1 12 


WR ITf (6.41) I 1 . Z ( 11 ) . X ( I M ) 

K 

1 13 


NXXXX=NP( IM ) 

K 

1 14 

38 

WRITF (6.42) ( (Y< IM. T ,L) . 1 = 1 . NV ) ,L = 1 .NXXXX) 

K 

1 15 

3Q 

RFTUPN 

K 

1 16 

r 


K 

1 17 

40 

FORMAT (///39H NO SOLUTION ON TAPE FOR THE COND I T 1 0NS//5X * 6H S/R= 

K 

1 18 


1F12.6.37X.9H EVALUATF/6X.2HZ ( I 1 . 2H ) = E1 6 . 8 , 38X . 3H Z< 11 ,1H)///) 

K 

119 

41 

FORMAT (1H12X.3H Z ( I 1 . 2H ) =F1 0.6/3X. 5H S/R=F 1 2 . 6// ! 20H T. DEG K 

K 

120 


1 LOG ( RHO ) LOGTP) LOG ( H/P ) GAMMA A/AA Z H 

K 

121 

2/RT Z* //) 

K 

122 

42 

FORMAT (2XF12.2.8F1 0.6 ) 

K 

123 


END 

K 

124- 


SUBROUTINE INTRP (N. X. Y. XINT.y INT ) 

L 

1 


DIMENSION X ( N ) « Y (N) 

L 

2 


Y T NT = 0 . 

L 

3 


DO 3 1=1 ,N 

L 

4 


SUMN= 1 . 

L 

5 


SUMD=1 . 

L 

6 


DO 2 J=1 « N 

L 

7 


IF ( J-I ) 1 424 1 

L 

8 

1 

SUMN=SUMN*(XINT-X(J ) ) 

L 

9 


SUMD=SUMD* (X ( I)-X ( J) ) 

L 

10 

2 

CONTINUE 

L 

1 1 

3 

Y I NT = Y I NT + Y ( I >*SUMN/SUMD 

L 

12 


return 

L 

13 


FND 

L 

14- 

— 4 

$ INP 

ITEST=1 . P5M= 1040,PT5M=1 . 746E+5 . U5M = 6 1 00 . RUN= 1 .RN=.0127$ 



$ INP 

ITEST=1 .P5M=1 040,PT5M= 1 . 746E + 5 . U5M = 6 1 00 , RUN= 1 . RN= . 0 1 27 , SAV = 0$ 



$ INP 

ITEST = 5.P5M= 1 04 0, RH0 C: M= .00485. U5M = 6 1 00 , RUN = 2 , RN= . 0 I 27* 



$ INP 

ITFST=5,P5M= 1 040 , RH05M= . 00485 , U5M=6 1 00 , RUN=2 , RN= . 0 1 27 , SAV=0* 



« INP 

ITEST = 3, P5M = 1040,PT5M= 1 . 746E+5 , QT5M = 2 . 743E + 7 , RUN = 3 , RN= . 0 1 27S 



$ INP 

ITEST = 3 4P5M = 1 040 4PTEM=1 . 746E + 5 * QT5M = 2 . 743E + 7 4 RUN = 3 4 RN = • 0 1 274SAV = 0$ 



— * 
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APPENDIX D — Continued 


A sample data printout is presented on the following page. The FORTRAN symbols 
in the printout and the corresponding symbols defined in the section entitled "Symbols" 
are as follows: 

PI 

Mm 

QTHOS 

q t (ref. 40) 

VI 

( V l)m 

QTFR 

q t (ref. 41) 

PT 

( p t)m 

QTPT 

k t (ref. 42) 

QT 

^t)m 

QTDE 

q t (ref. 43) 

RHOl 

( p l)m 

QTZO 

q t (ref. 44) 

RHOT 

( p t)m 

RN 

r g 

P 

P 



RHO 

P 



T 

T 




H h 

S/R s/R 

Z Z 

GAME y E 

A a 

V V 

M M 

NRE N Re 

QTCO q t (ref. 39) 


70 



RE AL-A I R DAT A R f. DUCT ION R°CGRaM OF MILLER 


ALL PHYSICAL QUANTITIES IN ”KS UNITS- NASA SP-7012 


MEASURED INPUTS 


RUN PI VI PT 

2.000F+00 1.040E+03 6.LCCR+03 0. 



R HOI RH3T TIME 

4.8S0E-03 0. 1 . 000E-04 


FREE-STREAM CONDI TIONS 

P RHO T H S/R Z GAME A V M NRE 

1.043E *0 3 4.850E-03 7«4fc2E*02 7.6355 *05 3.1775*01 1. 000 c *00 1,3615 *00 5.400E*02 6.1J 0F*0 3 1.130C *01 8 . 520E*05 


STATIC CONDITIONS BEHIND NORMAL SHOCK 

_P RHO T H_ S/R Z GAME _ A V M NRE 

1. 674E *05 6.225E-02 6 . 6 1 9 E + 0 3 L . 9 2 65 + 07 4.399E+0"l 1.414E+00 1.142E+00' 1.752E+03 4. 757E+02~ ~2. 7 1 5E-01 1.717E + 05 


STAGNATION CONDITIONS BEHIND NORMAL SHOCK 

p "pHO T H S/R ~l GAME A 

1, 746E+05 6.458F-02 6.643E+J3 1.9375+07 4.399E+01 1.416S+00 1.142E+00 1.757E+03 


STA GNATION POINT HEAT TRANSFER PRE DICTIO N S 

QTCO QTHOS QTFR OTPT QTDE OT ZO RN 

2.791E+07 2.345E+07 2.867F+C7 3.4475 + 07 2.893E+07 2.743F+07 1.270E-02 


APPENDIX D - Concluded 


APPENDIX E 


COMPUTER PROGRAM INPUTS 

The FORTRAN NAMELIST capability is used for data input with INP as the NAME- 
LIST name. The inputs necessary to utilize the computer program in appendix D are as 
follows: 


Program symbol 

Description and unit 

ITEST 

Identifies data reduction procedure to be used 

RUN 

Identifies facility test number 

T0LPT 

Desired tolerance of iteration involving p t 

T0LQT 

Desired tolerance of iteration involving q t 

T0LRH0 

Desired tolerance of iteration involving p^ 

SPH 

SPH = 0 for spherical q^. model 

SPH = 1 for transverse-cylinder q^. model 

Y0S 

Y0S = 0 for p t in q t predictions from Hansen (ref. 46) 
Y0S = 1 for p in q^. prediction from Yos (ref. 47) 

QD0 

QD0 = 0 if do not want various q^. calculated 
QD0 = 1 if want various calculated 

TWALL 

TWALL = 0 for T w = T amb 

TWALL = 1 if want to calculate T w using equation (11) 

TAU 

Test time, t, sec 

ETA 

Product PCpk of model surface material, W^-sec/m^-K^ 

TWE 

^w,amb> ^ 

RN 

r g, m 

P5M 

(Pl)m- N /“ 2 

PT5M 

(Pt)m> N/m2 
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Program symbol 

Description and unit 

U5M 

( V l)m’ m / sec 

QT5M 

(*t)m» W/m2 

RH05M 

(Pl) m > kg/™ 3 

RH0T5M 

(p t ) m , kg/m3 

SAV 

SAV = 0 for use of AEDC tape 

SAV = 1 for use of subroutine SAVE (K) 


In using a specific ITEST, only the measured inputs associated with that ITEST need 
be included. For convenience, the required inputs corresponding to a specific ITEST are 
indicated in the following table: 


ITEST 

P5M 

PT5M 

U5M 

QT5M 

RH05M 

RHOT5M 

1 

X 

X 

X 




2 


X 

X 

X 



3 

X 

X 


X 



4 

X 


X 

X 



5 

X 


X 


X 


6 

X 

X 



X 


7 


X 



X 

X 

8 

X 

X 




X 

9 


X 

X 



X 


To reduce the number of inputs in INP, values are assigned (within the computer 
program) to several input quantities; hence, these quantities may be excluded from INP, 
unless a change in the quantity from the assigned value is desired. The assigned values 
are as follows: 


73 



APPENDIX E - Concluded 


Program symbol 
T0LPT 
T0LQT 
T0LRH0 
SPH 
Y0S 
QD0 
TWALL 
TAU 
ETA 
TWE 
RN 
SAV 


Assigned value 
0.001 
0.005 
0.001 
0 
1 
1 
0 

1 x 10" 4 
2.045 x 10 8 
300 
0.01 
1 
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Figure 1.- Sketches illustrating expansion- tube cycle. Region (E) denotes 
free-stream conditions (i.e., test region). 
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Figure 2.- Comparison of several methods for predicting stagnation-point heat-transfer 
rate in air. r_. = 1.27 cm (sphere); T w = 300 K. 
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